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ABSTRACT 


A  series  of  experimental  studies  has  been  conducted  on 
transistor  elements  suitable  for  microcircuits.  A  relation¬ 
ship  of  transient  radiation  effects  to  active  volume  and  other 
parameters  has  been  made.  The  experimental  and  theoretical 
method  developed  offers  a  way  of  obtaining  minority  carrier 
lifetimes/  depletion  layer  width,  junction  area,  diffusion 
length,  and  active  volume  when  the  absorbed  dose  and  the  ion¬ 
ization  efficiency  are  known  in  the  transistor  element  of  a 
microcircuit.  Conversely,  given  the  above  parameters  it  is 
possible  to  determine  the  primary  photocurrent  of  the  micro- 
circuit  junction  under  a  variety  of  radiation  conditions.  De¬ 
velopments  have  been  made  on  the  Boltzmann  equation.  It  has 
been  shown  that  the  fundamental  Boltzmann  equation  describing 
charge  transport  in  semiconductors  can  be  used  to  modify  the 
usual  continuity  equation  to  obtain  a  more  accurate  circuit 
model  for  transistor  elements  in  microcircuits.  This  should 
result  in  changes  in  the  Ebers-Moll  model,  charge  control  model, 
and  the  Linvill  model  to  take  into  account  the  thin  layers  used 
in  microcircuit  diffusion  processes .  The  oxide-isolated  micro- 
circuits  have  been  tested  and  compared  with  the  usual  back- 
biased  isolation  junction  dxode  type.  Circuits  used  were  the 
same  except  for  the  type  of  isolation.  Several  mock-up  cir¬ 
cuits  were  used,  and  it  was  determined  that  the  mock-up  tech¬ 
nique  was  reasonably  accurate  for  radiation  effects  studies. 
Proper  elements  must  be  included  for  the  parasitic  elements 
that  are  inherent  in  microcircuits. 


iii 


CONTENTS 


Section 

I  TRANSIENT  X-RAY  RADIATION  EFFECTS  ON  p-n 

JUNCTIONS  SUITABLE  FOR  MICROCIRCUITS  . 

Introduction  . 

Initial  Interaction  of  X  Rays  with  Matter.  .  .  . 

Instrumentation  and  Equipment . 

Experimental  Results  . 

Conclusions . 

II  TRANSIENT  RADIATION  TESTING  OF  ELECTRONIC 

COMPONENTS  AND  INTEGRATED  CIRCUITS  . 

Introduction  . 

Motorola  MC-201  and  XC-201  Integrated  Circuits  . 
Fairchild  pL-903  Micrologic  Element  and 

Mock-up . 

TI  -  Series  51  Microcircuits  . 

Signetics  SU-315K . 

FXR  Test  of  Crysta Ionics  2N3084  FET . 

FXR  Test  of  Diffused  Resistor . 

Thin-Film  Nichrome  Resistors  . 

The  WSMR  Gamma  LINAC  Experimental  Results 

on  Hall  Devices . 

Summary  of  Conclusions  . 

Ill  PRELIMINARY  DISCUSSION  OF  ENGINEERING  MODELS 

FOR  THIN-LAYERED  DEVICES  ......  . 

• 

Introduction  to  Fundamentals  of  Diffusion 

Theory  .....  . 

The  Distribution  Function . 

The  Boltzmann  Transport  Equation  . 

Remarks . 

Introduction  to  Diffusion  Theory  and 

Semiconductor  Analysis  . 

The  Fundamental  Semiconductor  Equations . 

General  Transistor  Models . 

Subsection  Summary  . 

Introduction  to  the  Application  of 

Invariant  Imbedding . 

Time  Independent  Reflection  and  Transmission 

Coefficients  . 

Remarks  on  Time  Dependent  Invariant 

Imbedding . 

IV  GENERAL  CONCLUSIONS  AND  RECOMMENDATIONS . 


Page 


1 

1 

1 

13 

17 

29 


31 

31 

31 


36 

42 

33 

57 


61 

67 


69 


69 

70 
81 


82 

82 

86 

94 


95 

95 


102 

105 


APPENDIX  I .  109 

APPENDIX  II .  117 

APPENDIX  III .  123 

APPENDIX  IV .  129 

REFERENCES .  137 

DISTRIBUTION  .  l40 


v 


ILLUSTRATIONS 


rj— 


3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


18 

19 


Time  Dependence  of  Ionized  Electron 
Concentration  . 

Approximate  Minority  and  Majority  Carrier 
Densities  in  a  p-n  Diode  Before  and  Immediately 
After  Irradiation  . 

Simulation  of  Primary  Photocurrent  in  a 

p-n  Junction . 

Open  Base  Transistor  under  Irradiation . 

AFWL  600-KV  Flash  X-Ray  Facility . 

Schematic  Diagram  of  Transistor  Test  Circuit .  .  . 

Dose  as  a  Function  of  X-Ray  Tube  Voltage 
4  Inches  from  the  Front  of  the  Anode . 

Ip  as  a  Function  of  X-Ray  Tube  Voltage . 

Collector  Primary  Photocurrent  vs  Time  for  a 
2N  1613  . 

Radiation  Pulse  Shape  as  Determined  by 
Scintillator . 

Collector  Photocurrent  for  a  2Nl6l3,  Base 

Open  Configuration . 

MC-201  and  XC-201  Circuit  Diagram  . 

liL-903  Logic  and  Circuit  Diagrams . 

Mock-up  of  Fairchild  pL-903  . 

Transient  Response  of  the  ijL-903  "nand" 

Gate  and  Its  Mock-up . 

Circuit  Diagram  of  SN  512  and  SN  513 
Micrologic  Circuit . 

The  Response  of  SN  513  Microcircuit  and 
Mock-up  to  Flash  X-Ray  Pulses.  All  Time 
Scales  are  0.1  |isec/div.  The  Dose  was  500  mr. 
vcc  =  6  Volts.  The  50  f)  Cable  was  Terminated 
With  a  51  0  Resistor  in  Series  With  a  47  |j|if 
Capacitor  . 

Circuit  Representing  TI  SN  513  Made  of 
Conventional  Components  . 

Radiation  Pulses  of  SN  513  Mock-up  Made  of 
Conventional  Components  . 

Transistor  Circuit  Schematic . 

Cross  Section  of  Microcircuits . 

Equivalent  Circuit  Showing  Parasitic 
Transistors  . 


6 

10 

11 

14 

15 

16 
18 

*20 

22 

28 

32 

37 

38 

40 

44 


45 

46 

46 

48 

48 


20 

21 

22 


vi 


48 


ILLUSTRATIONS 
(con 't ) 


Figure  Pag 

23  Mock-up  Circuit  of  SN  5^3  Using  Conventional 
Components  and  p-n-p  Transistors  to  Represent 
Parasitic  Transistors.  The  Collector  Resistor 

is  Divided  into  Two  300  0  Resistors .  50 

24  A  Six-Resistor  and  Diode  Conventional  Component 
Mock-up  Circuit  of  TI  Microcircuit  SN  512 

(Only  Two  Inputs  Shown) .  51 

25  TI  SN  511  and  SN  511A  Schematic .  52 

26  SU-315K  Dual  "Nor"  Gate  Schematic .  54 

27  SU-315K  Test  Circuit .  54 

28  2N3084  Test  Circuit .  58 

29  Silicon  Diffused  Resistor  .  60 

30  Transient  Radiation  Response  of  a  Diffused 

Silicon  Resistor .  60 

31  Radiation  Effect  on  the  Ohio  Semiconductor 

HS-51  Hall  Device .  63 

32  Radiation  Effects  on  the  Ohio  Semiconductor 

HS-51  Hall  Device .  64 

33  Radiation  Effect  on  Beckman  335  Device .  65 

34  The  Fluxes  in  a  Stratified  Medium .  96 

35  Typical  Transient  Radiation  Response  of  a  p-n 


Diode  for  Various  Values  of  Lp,  Ln/  and  Wt •  •  .  .  116 

36  Characteristic  Results  of  Charges  Across  the 
Transition  Region  at  the  Base  Collector 


Junction  of  a  Grown  Junction  Transistor . 118 

37  X-Ray  Tube  and  Window  in  the  Screen  Room . 124 

38  Collimated  Radiation  Beam  of  AFWL  600-KV 

Flash  X-Ray  System . 125 

39  Energy  Spectrum  of  a  400  kvp  Flash  X-Ray  Pulse.  .  126 

40  Gamma  Flux  Equivalent  to  1  Roentgen -Hour  as 

a  Function  of  Gamma  Energy . 127 

41  Schematic  Diagram  of  Light  Emitting  Diode 

Test  Circuit . 130 

42  Shielding  Layout  for  Light  Emitting  Diode  Tests  .  131 

43  Peak  Transient  Current  in  GaAs  vs .  Bias 
Current  at  1-00  KV  on  Anode  of  X-Ray  Tube 

and  4"  from  Wall . 132 

44  Peak  Transient  Current  in  GaAs  vs.  X-Ray  Tube 
Voltage  for  Open  Circuit  with  Decoupling 

Capacitor . 134 


vii 


LIST  OF  TABLES 

Table  Page 

I  Physical  Parameters  for  Transistors  .  26 

II  LINAC  22  Mev  Electron  Result?  on  the  MC-201 

and  XC-201 . 35 

III  (iL-903  Electron  Radiation-Induced  Response . 4l 

IV  LINAC  Electron  Test  of  |iL-903  Mock-up . .  .  42 

V  SN  511  and  SN  511A  Radiation-Induced  Outputs.  .  .  53 

VI  SU-315K  FXR  Test  Results .  56 

VII  SU-315K  LINAC  Test  Results .  57 

VIII  Peak  Value  of  Gate  Current  from  2N3084 .  58 

IX  Peak  Value  of  Drain  Current  (ip) .  59 

X  Hall  Voltage  Change  (constant  Current)  During 

Irradiation  . .  62 

XI  Operating  Characteristics  of  the  600-KV 

Fexltron  Flash  X-ray  System .  123 


viii 


SECTION  I 


TRANSIENT  X-RAY  RADIATION  EFFECTS  ON  p-n  JUNCTIONS 
SUITABLE  FOR  MICROCIRCUITS 


1.  Introduction 

Basic  radiation  effects  phenomena  near  a  p-n  junction  are 
reviewed  with  the  desire  to  relate  them  to  microcircuit  design 
and  the  transient  radiation  environment.  Experiments  were 
performed  to  substantiate  theoretical  predictions  and  relating 
dose  to  primary  photocurrents /  secondary  photocurrents,  and 
active  volume.  Theoretical  calculations  which  are  applicable 
to  junctions  of  different  active  volumes  were  made  on  the 
time -dependent  primary  photocurrent.  The  theory,  together 
with  the  experimental  data,  offers  a  method  of  obtaining  mi¬ 
nority  carrier  lifetime,  depletion  layer  width,  junction  area, 
diffusion  length,  and  active  volume  when  the  absorbed  dose  and 
the  ionization  efficiency  are  known.  Conversely,  given  the 
above  parameters,  it  is  shown  that  it  is  possible  to  predict 
the  primary  photocurrent  of  a  p-n  junction  suitable  for  micro- 
circuit  design. 

2 .  Initial  Interaction  of  X  Rays  with  Matter 

a .  Ionization  of  Electrons 

When  X  rays  interact  with  matter  there  are  three  predomi¬ 
nant  mechanisms  of  interaction.  They  are  the  photoelectric, 
Comptoh,  and  pair  production  effects.  The  relative  importance 
of  each  of  these  depends  on  the  X-ray  energy  and  the  atomic 
number  of  the  absorbing  material. 

For  low  energies  (s:  100  kev)  and  large  atomic  numbers  the 
photoelectric  effect  is  the  predominant  process.  The  photo¬ 
electric  effect  is  the  process  whereby  the  energy  of  the  inci¬ 
dent  photon  is  totally  absorbed  in  ejecting  an  electron  from 
an  atom  and  imparting  kinetic  energy  to  it. 


1 


For  energies  between  100  kev  and  4  mev  the  Compton  process 
becomes  important .  It  represents  an  interaction  with  an  orbital 
electron  of  an  atom  to  produce  an  energy-degraded  photon  and  a 
recoil  electron. 

At  higher  energies,  above  1.02  mev,  the  pair  production 
process  becomes  increasingly  important  with  increase  in  photon 
energy.  It  is  an  absorption  process  for  X  rays  in  which  the 
incident  photon  is  annihilated  in  the  vicinity  of  the  nucleus 
of  an  atom,  and  an  electron-hole  pair  is  produced  (reference  l). 

When  an  X  ray  hits  the  semiconductor,  electron-hole  pairs 
are  produced  through  direct  ionization  by  photons  or  by  the 
slowing  down  of  electrons  and  positrons.  In  most  materials  an 
amount  of  energy  equal  to  two  to  four  times  the  ionization  po¬ 
tential  of  the  material  would  be  expended  per  ionization 
produced . 

In  gases  approximately  30  ev  are  dissipated  per  electron- 
ion  pair.  Gamma  rays  deposit  most  of  their  energy  via  secon¬ 
dary  electrons.  Hole-electron  production  is  directly  propor¬ 
tional  to  the  amount  of  energy  per  unit  volume  deposited  by 
both  primary  and  secondary  radiation. 

The  energy  spectrum  of  secondary  electrons  due  to  X-ray 
interaction  depends  upon  the  type  of  processes  involved  in 
their  generation  (photoelectric,  Compton,  or  pair  production). 
The  first  generation  of  electrons  will  have  predominantly  high 
energies.  The  next  generation  of  high-energy  secondary  elec¬ 
trons  results  again  in  a  l/(Energy)2  electron  spectrum.  If 
this  spectrum  still  has  some  high-energy  components,  these 
components  are  capable  of  producing  further  ionization.  These 
electrons  will  then  continue  to  lose  their  energy  by  further 
inelastic  scattering  until  finally  they  are  trapped  or  recom¬ 
bined  (reference  2). 

Ionization  within  a  semiconductor,  therefore,  results  in 
the  production  of  electron-hole  pairs.  Both  hole  and  electron 
pairs  may  carry  current  eithei  by  diffusion, when  there  is  a 
carrier  concentration  gradient,  or  by  drift  when  there  is  an 
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electric  field.  The  free  electron  will  be  attracted  towards 
the  positive  terminal,  and  the  hole  moves  towards  the  negative 
terminal  in  the  presence  of  an  applied  electric  field  (refer¬ 
ence  3 )  • 

If  the  air  surrounding  the  semiconductor  is  ionized,  the 
positive  ion  produced  can  move.  But  the  positive  ion  moves 
with  a  much  smaller  mobility  than  the  free  electron  of  the  air, 
and  it  does  not  add  much  to  the  flow  of  current.  The  electrons 
produced  by  ionization  of  gases  around  a  semiconductor  device 
often  contribute  significantly  to  <’  '  results  of  transient 
radiation  measurements  unless  preventive  steps  are  taken  (ref¬ 
erences  4  and  5 ) • 

b.  Immobilization  of  Electrons 

Sooner  or  later  the  excess  electron-hole  pairs  produced  by 
radiation  are  captured  and  can  no  longer  contribute  to  the 
transient  effects  of  the  circuit.  One  method  of  immobilization 
would  be  by  direct  recombination  of  the  electron  and  positive 
ions.  But  this  process  happens  very  rarely  as  the  probability 
of  an  electron  coming  within  range  of  attraction  of  the  posi¬ 
tive  ion  is  very  small.  Instead,  more  complicated  multistep 
processes  (trapping,  attachment,  and  defect  centers)  are  re¬ 
sponsible  for  electron  capture.  The  carriers  in  insulators 
can  combine  rapidly  (10-^  to  10”^  sec),  or  one  or  both  types 
of  carriers  may  become  trapped  for  a  longer  period  of  time  at 
a  crystal  defect  site . 

In  solids,  where  the  freed  electron  occupies  a  state  in 
the  conduction  band,  the  electron  returns  to  a  valence  band 
state  in  several  sequential  transitions  by  means  of  a  recom¬ 
bination  center.  These  recombination  centers  are  discrete 
levels  located  near  the  center  of  the  forbidden-energy  gap, 
and  are  so  called  because  they  aid  the  recombination  between 
electrons  and  holes.  This  process  simply  amounts  to  the  center 
first  capturing  an  electron  and  subsequently  capturing  a  hole. 

A  recombination  center  thus  has  a  large  capture  cross  section 
(probability  that  any  localized  level  can  capture  an  excess 
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carrier)  for  both  electrons  and  holes.  In  some  cases,  depend¬ 
ing  upon  the  atomic  origin  of  the  level,  the  level  has  a  high 
capture  cross  section  for  one  of  the  mobile  carriers,  but  an 
extremely  small  cross  section  for  its  mate.  This  level  is 
called  a  trap,  because  it  tends  to  trap  the  carrier  for  which 
it  has  the  larger  capture  cross  section  without  capturing  the 
other  kind  of  carrier.  Some  trapping  phenomena  have  been  ob¬ 
served  in  which  excess  mobile  carriers  have  been  held  for  hours, 
or  days,  as  in  the  formation  of  F  centers  in  alkali  halides 
(reference  4 ) . 


Figure  1.  Time  Dependence  of  Ionized  Electron  Concentration 

Figure  1  shows  the  characteristic  shapes  of  the  signal 
that  would  be  observed  for  ionizing-radiation  pulse  lengths 
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which  are  short  or  long  compared  with  the  electron  lifetime. 

If  the  electronic  lifetime  t  is  short  compared  with  the  dura¬ 
tion  and  rise  and  fall  of  the  radiation  pulse,  one  would  expect 
the  electron  concentration  n  to  follow  the  radiation  pulse  as 
in  figure  lb.  But  if  the  electron  lifetime  is  long  compared 
with  the  duration  of  the  radiation  pulse,  the  characteristic 
shape  should  be  as  shown  in  figure  lc. 

In  a  typical  germanium  or  silicon  transistor  the  recom¬ 
bination  time  is  usually  of  the  order  of  a  few  nanoseconds  to 
a  few  microseconds.  The  recombination  occurs  indirectly  by 
use  of  a  recombination  center.  The  rate  of  recombination  is  a 
function  of  the  concentration  of  these  recombination  centers. 

The  purer  the  material  the  longer  is  the  lifetime  of  excess 
carriers.  The  maximum  experimental  value  for  very  pure  single 
crystal  silicon  is  about  500  microseconds. 

c .  Processes  of  Current  Flow  in  a  Semiconductor 

There  are  two  different  kinds  of  mobile  carriers  present 
in  a  semiconductor  device  and  each  kind  can  carry  current  either 
by  drift  in  a  potential  gradient  or  by  diffusion  in  a  density 
gradient.  The  two  different  kinds  of  mobile  carriers  are  the 
negatively  charged  electrons  and  positively  charged  holes. 

The  drift  component  results  from  a  movement  of  carriers 
due  to  the  force  of  an  electric  field  on  a  charge.  Since  op¬ 
posite  charges  attract,  the  electrons  drift  towards  the  posi¬ 
tive  potential  and  the  holes  drift  towards  the  negative 
potential . 

The  diffusion  of  carriers  results  from  a  density  gradient , 
i.e.,  the  tendency  of  particles  to  move  from  a  region  of  high 
concentration  to  a  lower  concentration.  Intuitively,  one  can 
see  that  the  shorter  the  distance  in  which  the  electron  density 
changes  (i.e.,  the  greater  the  gradient  of  electron  density) 
the  greater  will  be  the  diffusion  current  (reference  6). 
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a.  Reverse  Biased 


b.  Forward  Biased 


Figure  2.  Approximate  Minority  and  Majority  Carrier  Densit 
in  a  p-n  Diode  Before  and  Immediately  After  Irradiation 
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d.  Effect  of  Ionizing  Radiation  on  a  p-n  Junction 

Irradiation  of  a  semiconductor  p-n  junction  produces  an 
excess  carrier  concentration.  With  constant  voltage  applied  to 
the  junction,  the  resulting  current  pulse  depends  on  the  radi¬ 
ation  pulse  width,  the  pulse  shape,  and  the  recombination  time 
of  the  minority  carriers.  When  the  recombination  time  is  long 
compared  to  the  radiation  pulse  width,  the  concentration  in¬ 
creases  at  a  rate  proportional  to  the  radiation  dose  rate  until 
the  pulse  is  shut  off.  Then  the  semiconductor  relaxes  to  its 
preirradiation  steady-state  value  within  some  characteristic 
time.  But  if  the  recombination  time  is  short  compared  to  rad¬ 
iation  pulse  width,  then  the  excess-carrier  concentration  is 
proportional  to  the  dose  rate  at  any  particular  time  and  fol¬ 
lows  the  shape  of  the  radiation  pulse. 

The  variation  of  majority  and  minority  carriers  in  the 
vicinity  of  a  one -dimensional  p-n  junction  is  shown  in  fig¬ 
ure  2.  Solid  lines  indicate  the  equilibrium  carrier  concen¬ 
tration  under  bias  before  irradiation,  and  dotted  lines  indi¬ 
cate  a  typical  carrier  concentration  immediately  after  irradi¬ 
ation  (reference  j).  In  a  reverse  bias  condition  the  minority 
carriers  constitute  the  reverse  leakage  current  (see  figure  2a). 
Because  the  holes  in  the  n-region  are  injected  into  the  p- 
region  and  electrons  are  injected  into  the  n-region,  the  direc¬ 
tion  of  current  flow  is  from  the  n-region  to  the  p-region. 

When  the  diode  is  irradiated  by  a  pulse  of  ionizing  radi¬ 
ation,  electron  hole  pairs  are  generated  uniformly  throughout 
the  device.  Since  there  is  a  built-in  electric  field  in  and 
near  the  junction  transition  region,  the  carriers  generated  in 
this  region  will  be  swept  across  the  junction  and  collected  in 
less  than  a  few  nanoseconds.  If  the  radiation  pulse  is  much 
wider  than  a  few  nanoseconds,  these  carriers  constitute  a  cur¬ 
rent  which  has  essentially  no  time  delay  relative  to  the  radi¬ 
ation  pulse.  For  this  reason  this  current  is  refer .ed  to  as 
the  prompt  photocurrent  component . 
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In  the  regions  outside  of  the  junction  transition  region 
there  is  a  transient  increase  in  the  minority  carrier  densities 
as  shown  by  the  dotted  lines  in  figure  2a.  As  indicated  by  the 
dotted  lines,  the  minority  carrier  density  gradients  at  the 
edges  of  the  junction  have  increased.  This  would  increase  the 
minority  carrier  diffusion  tendency  (i.e.,  density  gradient)  in 
the  vicinity  of  the  junction  and  more  carriers  would  diffuse 
towards  the  junction.  Most  of  the  excess  minority  carriers 
which  are  within  a  diffusion  length  of  the  junction  diffuse  to 
the  junction  and  contribute  to  a  transient  current.  On  the 
average,  excess  minority  carriers  greater  than  one  diffusion 
length  from  the  p-n  junction  generated  in  the  p-  and  n-regions 
do  not  reach  the  junction  before  they  recombine.  Therefore, 
they  do  not  contribute  significantly  to  the  transient  current. 

The  diffusion  component  of  current  generally  takes  a  time 
which  is  large  compared  to  the  duration  of  the  radiation  pulse 
to  diffuse  to  the  junction.  It  is  therefore  referred  to  as  a 
delayed  photocurrent  component.  The  excess  holes  in  the  n- 
region  diffuse  towards  the  p-region,  so  that  the  radiation- 
induced  delayed  component  of  current  is  from  the  n-  to  the  p- 
region.  Since  the  depletion  and  diffusion  photocurrent  is  in 
the  same  direction  as  conventional  current  flow  in  a  reverse 
biased  diode,  the  net  result  is  a  transient  photocurrent  super¬ 
imposed  on  the  steady-state  leakage  current. 

A  forward  biased  junction  can  be  regarded  as  a  short  cir¬ 
cuit  in  most  instances.  The  voltage  drop  due  to  applied  volt¬ 
age  across  the  junction  will  be  small  and  normally  will  not 
completely  cancel  out  the  built-in  junction  potential.  Thus 
an  electric  field  will  still  exist  in  the  same  direction  as  in 
a  reverse  biased  junction.  Carriers  which  are  generated  in 
the  junction  transition  region  will,  as  before,  be  swept  across 
the  junction  by  this  field  and  collected  within  a  few  nano¬ 
seconds.  This  current  will  be  in  opposition  to  the  normal  for¬ 
ward  bias  current .  As  indicated  by  the  dotted  lines  in  fig¬ 
ure  2,  the  majority  carrier  densities  immediately  after  irradi¬ 
ation  will  not  be  altered  appreciably,  but  the  minority  carrier 
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density  gradient  will  decrease.  The  net  effect  of  pulsed  X  rays 
on  a  forward  biased  junction  is  a  transient  photocurrent  oppo¬ 
site  in  direction  to  the  normal  steady-state  current. 

The  derivation  of  the  junction  transient  current  is  given 
in  appendix  I.  This  is  the  same  type  of  result  given  by  Brown, 
van  Lint,  Caldwell,  Keister,  et  al. (references  8,  9>  10,  and 
11).  The  photocurrent  derived  is  for  a  radiation  pulse  of  mag¬ 
nitude  g  and  duration  t  as  given  by  the  equation 


ipp(t)  =  qAg 


0  S  t  £  tQ  (la) 


i 


PP 


(t) 


(ib) 


In  appendix  II  it  is  shown  that  the  width  of  the  depletion 
region  is  a  function  of  the  bias  voltage  (reference  12).  Thus 
the  primary  photocurrent  consists  of  the  following  two  compo¬ 
nents:  (l)  A  voltage  insensitive  diffusion  component  that  is 

contributed  to  by  the  p  and  n  regions;  and  (2)  a  voltage  sen¬ 
sitive  drift  component  that  is  contributed  to  by  the  depletion 
layer.  Usually  the  voltage  insensitive  diffusion  component  is 
much  greater  than  the  voltage  sensitive  drift  component  (ref¬ 
erence  13)*  The  primary  photocurrent  can  be  considered  as  a 
constant -current  generator  connected  across  the  junction  as 
shown  in  figure  3* 

If  the  collector-base  region  of  a  planar,  discrete,  or 
integrated  circuit  transistor  is  used  as  a  diode,  then  most  of 
the  primary  photocurrent  usually  originates  in  the  collector 
region.  With  the  above  approximations  of  neglecting  the  de¬ 
pletion  region  and  the  base  diffusion  component,  the  transient 
photocurrent  response  to  a  radiation  pulse  is  of  the  form 


where 

.  -IQ 

q  =  electron  charge  =  1.6  x  10  coulombs 

Ac  =  collector -to-base  junction  area 

g  =  electron -hole  pair  generation/cm^-sec 

Lc  =  collector  minority  carriers  diffusion  length 

tc  “  lifetime  of  minority  carriers  in  the  collector 

t  =  X-ray  pulse  width 

u(t)  =  unit  step  function  starting  at  t  =  0 
u(t  -  tQ)  =  unit  step  function  starting  at  t  =  tQ 


Figure  3.  Simulation  of  Primary  Photocurrent 

in  a  p-n  Junction 


The  generation  rate,  g,  for  a  silicon  transistor  can  be 
calculated  from  (references  14  and  15 ) 

y  =  4.2  x  1015  y  pairs/cm^-sec  (5) 


9  =  — 


(E  1.6  x  10* 


10 


fe=. t 


density  of  silicon  =  2.42  g/cm^ 
average  energy  in  electron  volts  to  produce  an 
electron  hole  pair  «  3*6  ev  for  silicon 
dose  rate  in  rads/sec  (Si) 

Effect  of  Ionizing  Radiation  on  a  Transistor 

Consider  an  n-p-n  transistor  connected  in  the  open  base 
configuration  and  biased  as  shown  in  figure  4.  When  the  tran¬ 
sistor  is  irradiated  by  a  short  pulse  of  ionizing  radiation, 
electron-hole  pairs  are  created  uniformly  through  the  volume. 
Just  as  before  in  the  case  of  the  diode,  the  radiation  gener¬ 
ated  electrons  and  holes  will  drift  across  the  junction  and 
diffuse  to  the  junction.  This  will  result  again  in  a  prompt 
photocurrent  component  from  the  transition  region  and  a  delayed 
component  from  the  diffusion  of  minority  carriers  from  one 
region  to  another.  With  an  n-p-n  transistor  the  electrons  in 
the  base  will  diffuse  to  the  collector-base  junction,  and  they 
will  be  collected  as  a  component  of  the  delayed  photocurrent. 


where 

o  = 

E  = 

Y  = 
e. 
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Figure  4.  Open  Base  Transistor  Under  Irradiation 


In  the  absence  of  these  electrons,  the  holes  left  behind 
form  an  excess  positive  charge  in  the  base.  At  the  time  the 
electrons  are  diffusing  from  the  base,  holes  are  diffusing 
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from  the  collector  to  the  collector-base  junction,  and  they  are 
swept  into  the  base  region.  These  holes  will  further  increase 
the  excess  positive  charge  in  the  base.  The  excess  majority 
carriers  (holes)  in  the  base  will  increase  the  forward  bias  of 
the  emitter-base  junction,  which  causes  more  electrons  from  the 
emitter  to  flow  into  the  base.  The  electron  current  entering 
the  base  from  the  emitter  is  determined  by  the  current  gain  of 
the  device  as  well  as  the  radiation-induced  charge  in  the  base. 
It  results  in  a  collector  current  transient  (secondary  photo¬ 
current)  in  addition  to  the  collector  component  of  the  primary 
photocurrent.  These  injected  minority  carriers  (electrons  from 
the  emitter  into  the  base)  will  continue  to  diffuse  to  the 
collector-base  junction  and  drift  across  the  junction  as  long 
as  an  excess  positive  charge  remains  stored  in  the  base  region 
(references  7  and  10). 

From  the  above  it  is  seen  that  the  photocurrent  consists 
of  the  following  contributions:  (l)  A  voltage-sensitive  deple¬ 
tion  layer  component  contributed  by  the  emitter -base  and 
co  Hector -base  junction,  and  (2)  a  voltage-insensitive  diffu¬ 
sion  component  contributed  by  the  emitter,  base,  and  collector 
regions . 

In  planar  and  mesa  transistors  most  of  the  primary  photo¬ 
current  usually  originates  in  the  collector  region  (reference 
16) .  The  collector  diffusion  length,  L  ,  is  much  greater  than 
the  base  width  and  the  depletion  width.  This  means  that  the 
number  of  holes  available  from  the  diffusion  of  holes  from  the 
collector  to  base  is  much  greater  than  the  diffusion  of  elec¬ 
trons  from  the  base  to  collector.  Thus,  after  the  termination 
of  the  radiation  pulse,  holes  generated  in  the  collector  region 
within  one  diffusion  length  of  the  junction  will  diffuse  into 
the  base.  The  charge  accumulated  in  the  base  due  to  diffusion 
of  carriers  from  the  collector  will  therefore  not  reach  a  max¬ 
imum  until  approximately  a  time  equal  to  one  to  two  times  the 
lifetime  of  minority  carrier  in  the  collector  (t  ). 
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It  is  postulated  that  the  substrate  collector  junction  will 
contribute  only  primary  photocurrent  to  the  total  photocurrent. 
The  substrate  to  base  dimension  is  wider  than  a  diffusion  length 
so  that  the  substrate-collector-base  transistor  has  a  very  low 
current  gain  and  poor  transistor  action.  It  will  then  follow 
that  the  secondary  photocurrent  will  be  small. 

3 .  Instrumentation  and  Equipment 

a .  General  Layout 

Figure  5  shows  the  layout  of  the  flash  X-ray  facilities  at 
Kirtland  Air  Force  Base.  The  controls  for  the  flash  X-ray  sys¬ 
tem  are  situated  outside  of  the  screen  room.  The  test  equip¬ 
ment  and  the  test  sample  are  enclosed  in  a  double-walled  rf 
shield  room.  Inside  the  room  the  rf  noise  level  is  reduced  by 
120  db.  See  appendix  III  for  a  more  detailed  discussion  of  the 
flash  X-ray  system. 

b .  Test  Circuit 

The  test  circuit  is  shown  in  figure  6.  Preliminary  test¬ 
ing  showed  that  in  order  to  keep  the  circuit  time  constant 
small  compared  to  the  transistor  response,  the  viewing  resistor 
had  to  be  less  than  200  ohms .  It  was  decided  to  use  a  viewing 
resistor  which  would  properly  terminate  the  coaxial  cable  with 
its  characteristic  impedance  of  50  ohms . 

In  order  to  minimize  the  cable  length,  and  hence  the  con¬ 
stant  of  the  viewing  circuit,  an  amplifier  was  used.  This  al¬ 
lowed  the  leads  to  be  shortened  to  approximately  2  feet.  The 
amplifier  and  power  supply  were  shielded  from  incident  radia¬ 
tion  by  a  wall  made  of  2-inch  thick  lead  bricks. 

Air  ionization  will  cause  a  shunt  leakage  path  on  the 
order  of  10  megohms  (reference  5)*  Since  the  viewing  resistor 
is  very  much  smaller  than  the  shunt  resistance  due  to  air 
ionization,  the  effect  of  air  ionization  was  neglected.  In 
order  to  eliminate  any  effect  of  radiation  on  the  cable  or 
resistors,  a  1-inch  thick  lead  block  was  placed  between  the 
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X-ray  source  and  the  circuit.  A  1-inch  hole  was  drilled  in  the 
block  to  allow  the  X  rays  to  impinge  upon  the  target  transistor. 

c .  Noise  Level 

The  ambient  noise  level  is  due  primarily  to  the  noise  gen¬ 
erated  in  the  amplifier.  The  amplifier  noise  level  is  specified 
as  less  than  50  microvolts  peak  to  peak. 

The  transient  noise  level  is  caused  by  the  operation  of  the 
flash  X-ray  system.  A  second  screen  room  would  probably  reduce 
the  transient  noise  level  sufficiently  that  it  would  be  masked 
by  the  ambient  noise  from  other  sources. 

Noise  level  is  particularly  important  when  irradiating 
devices  with  small  active  volumes  and  fast  response  times. 

Small  active  volumes  will  give  low  amplitude  response,  and  fast 
response  times  will  force  small  viewing  resistors  to  keep  the 
time  constant  down.  The  combination  of  these  two  factors 
forces  use  of  maximum  gain  of  the  amplifier  and  oscilloscope 
preamplifiers.  The  use  of  maximum  gain  of  the  amplifier  is  the 
worst  condition  from  a  noise  standpoint.  During  recent  tests, 
the  response  of  a  2N917  transistor  was  masked  by  noise  so  that 
only  crude  order-of-magnitude  results  could  be  obtained. 

4 .  Experimental  Results 

a .  Dose  Dependence 

Tests  were  run  using  the  circuit  shown  in  figure  6  to 
determine  peak  collector  current  as  a  function  of  dose  in 
milliroent gens  .  The  variation  in  dose  was  obtained  by  vatying 
the  anode  voltage  on  the  flash  X-ray  tube. 

A  plot  of  dose  vs .  X-ray  tube  voltage  is  shown  in 
figure  7.  The  resultant  curve  was  linear  with  a  constant 
slope  in  the  region  of  interest.  The  curve  can  be  approxi¬ 
mated  by  equation  4. 

mr  =  0.572V  +  29  (4) 
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where 

mr  =  dose  in  milliroentgens 
V  =  X-ray  tube  voltage  in  kilovolts 

As  a  result  of  this  experiment  all  further  tests  were  run  using 
anode  voltage  as  an  approximate  measure  of  dose. 

A  plot  of  peak  primary  photocurrent  vs.  anode  voltage  for 
a  2Nl6l3  transistor  is  shown  in  figure  8.  The  slight  flatten¬ 
ing  of  the  curve  above  450  KV  can  be  attributed  to  two  factors: 

(1)  The  transistor  response  is  to  some  extent  sensitive  to  dose 
rate.  As  the  anode  voltage  is  increased  the  shape  of  the  X-ray 
pulse  is  distorted.  This  distortion  is  much  more  pronounced  at 
high  anode  voltages  and  causes  a  change  in  dose  rate  which  is 
less  than  the  expected  linear  change  with  anode  voltage. 

(2)  The  transistor  may  be  approaching  a  hole-electron  pair  pro¬ 
duction  rate  saturation. 

As  a  result  of  this  experiment,  testing  for  variation  with 
other  parameters  was  done  in  the  linear  portion  of  the  response 
curve;  400  KV  anode  voltage  was  selected  to  give  maximum  dose 
and  still  remain  in  the  desired  portion  of  the  curve. 

b .  Active  Volume 

The  active  volume  depends  upon  junction  area,  diffusion 
length,  and  depletion  width.  All  of  these  are  constant  except 
depletion  width,  which  is  a  function  of  applied  voltage.  An 
experiment  was  set  up  to  determine  the  effect  of  bias  voltage 
upon  the  peak  collector  current . 

Using  the  circuit  shown  in  figure  6,  the  voltage  was 
varied  in  steps  from  0  to  40  V  DC.  A  plot  of  log  Ip  vs.  log 
bias  voltage  resulted  in  a  straight  line,  where  Ip  is  the  peak 
transient  current  due  to  an  X-ray  pulse.  From  the  slope  of  the 
line  the  relationship  between  Ip  and  bias  voltage  is  given  by 


KV1/18  +  I 


(5) 
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Figure  9*  Collector  Primary  Photocurrent  vs  Time  for  a  2N  1613 
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where 

I  =  peak  collector  current  in  ma 
P 

V  =  bias  voltage  in  volts 

K  =  constant 

I  =  zero  bias  current 
o 

This  result  shows  an  apparent  conflict  with  the  theoretical 
results  derived  in  appendix  II.  The  conflict  can  be  resolved  by 
considering  the  relative  times  and  magnitudes  of  the  two  com¬ 
ponents  of  the  primary  photocurrent . 

The  drift  component  of  the  primary  photocurrent  is  de¬ 
veloped  in  the  depletion  region.  The  diffusion  component  comes 
from  the  region  up  to  one  diffusion  length  beyond  the  depletion 
region.  In  a  semiconductor  with  a  moderately  long  minority 
carrier  lifetime,  a  diffusion  length  is  much  longer  than  the 
depletion  region  width.  Therefore,  the  contribution  from  the 
diffusion  current  is  much  larger  than  that  from  the  drift 
current . 

The  diffusion  current  is  of  the  form  erf  ,/t/r .  The 

time  to  peak  current  is  a  function  of  t,  the  minority  carrier 
lifetime.  The  drift  current  is  of  the  form  I^gCt)  where  g(t) 
is  the  generation  rate.  The  drift  current  will  follow  the 
radiation  pulse  with  no  delay. 

Figure  9  shows  theoretical  and  experimental  curves  of 
photocurrent  plotted  against  time.  When  the  theoretical  curve 
for  diffusion  current  is  subtracted  from  the  experimental 
curve,  the  result  is  the  drift  component  of  the  photocurrent. 

The  drift  component  has  followed  the  radiation  pulse  as  can  be 
seen  by  a  comparison  with  figure  10  which  shows  the  shape  of 
the  radiation  pulse  as  determined  by  a  scintillator. 

The  peak  of  the  drift  component  does  not  coincide  with  the 
peak  of  the  diffusion  component  of  the  photocurrent.  The  peak 
of  the  sum  lies  somewhere  in  between .  Measurement  of  the  peak 
photocurrent  is,  therefore,  not  an  accurate  measurement  of  the 
change  in  amplitude  of  the  drift  component  of  the  photocurrent. 
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The  drift  components  of  the  photocurrent  for  two  different 
bias  levels  were  separated  from  the  total  photocurrent.  Their 
average  values  were  determined  to  be  approximately  70  ua  for  a 
bias  of  10  volts  and  45  |ia  for  a  bias  of  2  volts.  When  the 
points  were  plotted  on  log  log  graph  paper,  the  slope  of  the 
line  was  found  to  be  0.275*  This  result  may  be  stated  as 

1/3.64 
d 

where  1^  is  the  drift  component  of  the  photocurrent  which  com¬ 
pares  favorably  with  the  result  from  appendix  II,  equation  ( 162 ) 

ld  =  R  v1/?  (7) 


(6) 


c .  Transistor  Response 


Two  transistor  types  were  experimentally  investigated. 

They  were  the  2Nl6l3  and  2N91T*  The  transistors  were  operated 
with  a  10-volt  collector  bias  voltage.  Data  were  obtained  for 
a  base  emitter  short  circuit  configuration  as  shown  in  figure  6 
and  for  a  base  open  circuit  configuration. 

The  intent  was  to  obtain  comparative  results  for  the  two 
transistor  types.  This  intent  was  not  realized  because  of  the 
noise  problems  discussed  earlier. 

( 1 )  Base  Emitter  Short 

Equation  (la)  will  give  an  approximate  theoretical  expres¬ 
sion  for  the  primary  photocurrent  if  the  base  component  is 
neglected,  as 


qAc  g[Wt  +  Lp  erf  0  s  *  s  fco 


(8) 


The  prompt  current,  ip(t),  can  be  separated 
component,  i^(t);  they  are  given  from  above 


from  the  diffusion 
as 
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ip(t )  =  q  Ac  g  Wt 


0  <  t  <  t 


(9) 


i,(t)  =  *3  A  g  L  erf  *Jt/r  0  t  £  t 


(10) 


In  figure  9  the  experimental  depletion  component  (prompt  or 
drift  component)  was  separated  from  the  diffusion  component 
(delayed  component).  The  diffusion  component  data  can  be 
fitted  with  an  expression  given  by 

i  (t)  =  q  A  g  L  erf  *Jt/j  =  (400  Ma)  erf  Jt/ .J>  psec 


0  £  t  £  t 


(id 


for  a  type  2Nl6l3  transistor  col lector -base  junction. 

Since  L  =  ■Jfyr ,  one  can  calculate  out  the  value  of  L  by 

r 

choosing  a  typical  value  for  the  diffusion  constant 
(Dp  =  13 .3  cm^/sec) 


T 


LP  =  '/Vp=  ( .3  x  10~6  sec ) 


(12) 


Lp  as  2  x  10-5  cm 


(13) 


From  equation  (10)  and  figure  9  it  is  seen  that 


q  Ac  g  Lp  =  400  na 


(14) 


Or,  solving  for  q  Ac  g,  one  may  obtain 

(13 

Using  the  value  calculated  for  Lp  in  equation  ( 13 )  equation  (l4) 
becomes 


q  Ac  g  =  ^ 
P 
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HUM  -  .  -i — ■ 


q  a 


„  =  400  Ma 

c  2  x  10'5 


=  200  — 
cm 


(16) 


The  shape  of  the  experimental  depletion  current  followed 
the  shape  of  the  X-ray  pulse.  This  pointed  to  the  fact  that  the 
radiation  pulse  was  not  constant  for  the  0.2  Msec.  The  varia¬ 
tion  in  the  radiation  rate  can  be  compensated  for  by  taking  the 
average  value  of  the  depletion  component  as  an  approximation. 
Using  an  average  value  of  about  70  ua  for  an  applied  reverse 
bias  of  10  v  across  the  diode  composed  of  the  collector-base 
junction  of  the  2Nl6l3#  one  can  solve  for  an  average  value  of 
the  depletion  width  .  Solving  for  out  of  equation  (9)  one 
obtains 


uo 


W  =  — C - 

t  q  Ac  g 


(17) 


Using  the  results  of  equation  (l6)  and  an  average  value  for 
ip(t),  the  average  value  for  Wfc  becomes 


wt  = 


'0  ua 


raa/cm 


=  3-5  x  10' 


(18) 


This  value  of  the  depletion  width  is  about  the  same  value  one 
would  theoretically  calculate. 

If  an  average  energy  of  3*6  ev  is  required  to  create  an 

electron-hole  pair  in  silicon,  a  dose  rate  of  1.4  x  10^  r/sec 

19 

will  produce  a  generation  rate  of  about  7  x  10  electron-hole 
pairs/cm^-sec .  Equation  (l6)  gave 


q  A  g  =  200  §  =  .2  as* 

c  cm  cm 


(19) 


If  q  =  1.6  x  10” ^  coulombs  and  g  =  7  x  lO1^,  one  can 


solve  for  A  . 

c 
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&  _  0 . 2  _  _ 0 . 2 _ 

c  qg  (1.6  x  10-1^)(7  x  1019) 

Ac  =  0.01785  =  17.85  x  10"5  cm2  (21) 

In  table  I ,  the  values  which  were  able  to  be  predicted 
from  radiation  measurements  and  the  measured  transistor  values 
are  tabulated. 


Table  I 

Physical  Parameters  for  Transistors 


Measured 

2N1613 

Radiation 

Predicted 

2N1613 

Measured 

2N917 

Collector  Impurity 
Concentration 

5x10^  cm”^ 

2x10^  cm" 

Base  Impurity 
Concent  rat  ion 

101?  cm-5 

5x10^  cm” 

Emitter  Impurity 
Concentration 

1020 

1020 

Collector  Base 
Junction  Area 

5x10”^  cm2 

17.85xl0"5  cm2 

5x10”^  cm2 

Emitter  Base 
Junction  Area 

-4  2 

5x  10  cnr 

6x10”^  cm2 

Collector  Lifetime 

0.5  lisec 

0.3  (asec 

Base  Width 

2  microns 

0.5  micron 

Chip  Thickness 

{Xy)  jei  crons 

80  microns 

Co  Hect  or -Base 
Depletic?'  Width 
at  volts 

•2.5 

Microns 

5.5 

microns 

The  discrepancy  the  radiation  predicted  and  the 

measured  value"  irefiecr  uncertainties  in  measuring  the  ab- 

sorbed  dose,  measuring  ahysical  parameters  of  the  transistor, 

experiments*  *rro»  ""d  the  approximations  used  in  the  analysis, 
such  as  the  y  T  the  junction. 
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By  looking  at  table  I  again,  it  can  be  seen  that  the 
collector -base  junction  area  of  a  2Nl6l3  is  a  hundred  times 
larger  than  the  2N917.  If  all  the  other  parameters  remained  the 
same,  a  reduction  of  a  hundred  would  be  theoretically  predicted 
from  equation  (l)  in  the  peak  transient  current  of  the  2N917 
over  the  peak  transient  current  of  the  2N1613.  Experimental  re¬ 
sults  showed  a  reduction  of  approximately  forty. 

(2 )  Base  Open 

Figure  11  shows  the  collector  current  wave  form  of  a  2Nl6l3 
transistor  that  was  irradiated  in  the  base  open  configuration. 
The  current  consists  of  primary  and  secondary  photocurrent.  The 
amplitude  of  the  primary  photocurrent  is  approximately  the  same 
as  the  amplitude  of  the  photocurrent  obtained  in  the  base 
shorted  configuration. 

In  a  transistor  whose  drift  component  is  small  compared  to 
the  diffusion  component  of  the  primary  photocurrent,  the  magni¬ 
tude  of  the  secondary  photocurrent  can  be  approximated  by 


I 


sp 


(22) 


where 

Isp  -  peak  secondary  photocurrent 

hfe  =  small  signal  current  gain  measured  at  I 

ipp  =  peak  primary  photocurrent . 

This  approximation  is  good  for  open  base  configuration  only  and 
does  not  contain  any  time  information. 

The  secondary  photocurrent  is  attributed  to  a  positive 
charge  build-up  in  the  base.  It  remains  until  the  charge  in 
the  base  is  dissipated.  An  experiment  was  run  to  determine 
what  impedance  between  base  and  emitter  is  necessary  to  elimi¬ 
nate  the  secondary  photocurrent .  Starting  with  the  base  emitter 
impedance  of  zero,  the  impedance  was  increased  to  100  kilohms 
without  observing  any  secondary  photocurrent  for  doses  up  to 
250  mr.  One  hundred  kilohms  is  much  larger  than  any  resistors 
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used  in  molecular  integrated  circuits.  This  would  imply  that 
most  transistors  used  in  integrated  circuits  may  be  treated  as 
if  their  base  emitters  were  shorted. 

5 .  Conclusions 

The  following  conclusions  can  be  made  from  the  theoretical 
and  experimental  studies  that  were  carried  out:  (l)  The  active 
region  should  be  kept  as  small  as  possible.  This  implies  mini¬ 
mization  of  emitter -base  junction  and  collector  junction  areas, 
minority  carrier  lifetimes,  and  depletion  widths.  The  base  col 
lector  junction  area  can  be  minimized  by  using  epitaxial  tech¬ 
niques  and  selecting  the  smallest  possible  device  consistent 
with  the  required  power  or  current  levels.  Minority  carrier 
lifetime  can  be  reduced  by  gold  doping  and  using  high-frequency 
devices.  The  depletion  width  can  be  minimized  by  using  a  mini¬ 
mum  collector  bias  voltage.  (2)  The  circuit  should  use  the 
lowest  impedances  consistent  with  design  criteria.  The  current 
source  behavior  of  the  radiation  effect  indicates  the  desirabil 
ity  of  using  low  impedances.  The  secondary  photocurrent  which 
is  a  result  of  a  high  base  emitter  impedance  is  another  factor 
in  the  desire  for  low  impedances.  (5)  The  experimental  data 
with  the  theory  developed  here  offer  a  method  of  obtaining 
minority  carrier  lifetimes,  depletion  layer  width,  junction 
area,  diffusion  length,  and  active  volume  when  the  absorbed 
dose  and  the  ionization  efficiency  are  known.  Conversely, 
given  the  above  parameters,  it  is  possible  to  determine  the 
primary  photocurrent  from  the  junction.  Thus,  it  is  possible 
using  these  techniques  to  make  some  general  recommendations 
for  design  and  fabrication  of  integrated  circuits  for  use  in  a 
transient  radiation  environment. 

In  current  molecular  integrated  circuit  design,  many  of 
these  features  that  are  desired  for  radiation  hardening  are 
being  incorporated  for  other  reasons .  Minimizing  transistor 
size  will  increase  manufacturing  yield.  Minimizing  lifetime 
leads  to  higher  frequency  devices.  Gold  doping  tends  to  reduce 
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inverse  a  and  therefore  leakage  current  in  some  transistor  con¬ 
figurations.  Small  valued  resistors  take  up  less  space  and  thus 
tend  to  increase  manufacturing  yield. 
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SECTION  II 


TRANSIENT  RADIATION  TESTING  OF  ELECTRONIC 
COMPONENTS  AND  INTEGRATED  CIRCUITS 

1.  Introduction 

Concurrent  with  the  theoretical  and  experimental  studies 
described  in  the  previous  section,  rad-  cion  tests  on  numerous 
components  and  integrated  circuits  were  conducted  at  the  Air 
Force  Weapons  Laboratory  flash  X-ray  facility  and  on  the  WSMR 
LINAC.  The  purposes  of  these  tests  were  (1)  to  attempt  to 
isolate  factors  in  the  design  and  fabrication  of  integrated 
circuits  which  have  the  greatest  effect  on  the  circuit  response 
during  a  pulse  of  ionizing  radiation;  and  (2)  to  evaluate 
several  special  thin-film  devices  which  appeared  to  be  less 
affected  by  the  radiation  than  integrated  circuit  active 
devices.  The  results  of  these  tests  are  presented  in  this 
section. 

2 .  Motorola  MC-201  and  XC-201  Integrated  Circuits 
a.  Background  Information 

The  Motorola  MC-201  and  XC-201  are  monolithic  silicon, 
epitaxial,  passivated,  integrated  circuits.  They  are  used  as 
four  input  "nand"  or  "nor"  gate  circuits  for  high-speed,  low- 
power  computer  systems.  Their  attractive  features  include 
high  stray  noise  immunity,  fast  switching  speeds,  and  low -power 
dissipation.  The  XC-201  and  the  MC-201  are  quite  similar  in 
operation  except  that  the  XC-201  has  a  much  faster  switching 
time  response.  The  MC-201's  have  switching  time  turn-on  delays 
and  turn-off  delays  of  60  and  50  nsec,  while  the  XC-201's  are 
respectively  8  and  17  nsec.  The  MC-201  and  XC-201  average 
propagation  delays  are  respectively  50  and  17  nsec.  The 
increased  speed  of  the  XC-201  is  achieved  by  the  use  of  the 
EPIC  process.  In  the  usual  monolithic  semiconductor  circuits, 
isolation  between  various  elements  on  the  same  substrate  is 
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accomplished  by  means  of  a  diffusion  which  yields  back-to-back 
p-n  junctions  between  different  elements.  With  the  application 
of  bias  voltage,  these  junctions  represent  reverse  biased  diodes 
with  a  very  high  back  resistance  which  provides  adequate  DC 
isolation.  Unfortunately,  as  each  p-n  junction  is  also  a 
capacitance,  there  remains  capacitive  coupling  between  compo¬ 
nents  and  the  substrate  at  high  frequencies.  The  EPIC  process 
completely  isolates  all  components  from  each  other  without  the 
use  of  back -biased  isolation  junction  diodes.  Isolation,  both 
electrical  and  physical,  is  achieved  by  a  layer  of  nonconductive 
material  which  completely  surrounds  and  separates  the  components 
from  each  other  and  from  the  common  substrate.  This  insulating 
layer  can  be  made  thick  enough  that  its  associated  capacitance 
is  negligible. 

From  the  transient  radiation  sensitivity  standpoint  the 
EPIC  process  offers  promise  of  increased  hardening  because  of 
the  elimination  of  the  highly  radiation-sensitive  back-biased 
diodes  for  isolation.  The  main  purpose  of  the  MC-201  and  XC-201 
transient  radiation  tests  was  to  determine  what  gains  (if  any) 
in  radiation  hardening  have  been  achieved  over  the  MC-201  by 
the  use  of  the  EPIC  process  in  the  XC-201. 


In- figure  12  are  shown  the  logic  and  approximate  circuit 
diagrams  for  the  MC-201  and  XC-201. 


v3,4,6  =  +  8  volts 

=  -  8  volts 

Figure  12.  MC-201  and  XC-201  Circuit  Diagram 
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The  operation  of  the  circuit  as  a  "nand"  gate,  for  example, 
is  as  follows:  When  any  of  the  nputs  are  at  ground  potential, 
the  voltage  at  point  "b"  is  also  at  ground  potential,  except 
for  the  small  positive  diode  forward  voltage  drop.  The  voltage 
at  point  "c"  is  two  diode  voltage  drops  below  point  "b"  or 
negative.  This  reverse  biases  the  transistor's  emitter-base 
junction,  and  thus  the  gate  is  turned  off.  When  all  of  the 
inputs  have  positive  applied  voltages,  the  gate  is  turned  on 
and  the  potentials  at  points  "b"  and  "c"  increase  to  provide 
base  current  to  the  transistor.  This  drives  the  transistor  into 
saturation  and  thus  reduces  the  voltage  at  "d"  to  a  small  value. 

.  Vi* 

b.  FXR  Tests  of  MC-201  and  XC-201 

The  MC-201  and  XC-201  were  tested  at  the  AFWL  flash  X-ray 
system  at  Kirtland  AFB,  New  Mexico.  The  tests  were  performed 
at  the  400-KV  level  with  a  0.2-nsec  pulse  width.  This  yields  an 
approximate  dose  rate  of  2  x  10^  roentgens/sec. 

In  figure  12,  the  external  test  circuit  and  conditions 
were  as  follows: 

a)  a  510-0  resistor  was  connected  from  terminals  6  to  5 

b)  Vg  =  +4  volts 

c)  Vj.  =  volts 

d)  Inputs  7#  9#  and  10  open 

e)  v0  =  0  or  +4  volts 

The  XC-201  was  tested  during  irradiation  with  V  both  zero 

o 

and  plus  four  volts.  No  transient  radiation  response  was  detect¬ 
able  above  the  noise  level  when  using  an  oscilloscope  gain  of 
0.005  volt/cm  with  a  peak-to-peak  noise  voltage  present  of 
approximately  0.0015  volt.  It  is  concluded  that  the  XC-201  is 
highly  radiation  resistant  to  gamma  doses  up  to  at  least 
2  x  106  roentgens/sec. 
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The  testing  of  the  MC-201  (under  identical  test  conditions 
as  for  the  XC-201)  did  yield  a  transient  radiation  response 
output  (tyj) .  These  transient  outputs  were  of  approximately 
the  same  duration  as  that  of  the  gamma  radiation  pulse,  and 
the  peak  values  were: 

a)  V^(peak)  *  -0.0035  volt  when  V0  =  0  volts, 

and 

b)  V^ipeak)  -  +0.0025  volt  when  V0  =  +4  volts. 

The  difference  in  polarity  for  the  two  states  is  most 
likely  attributable  to  the  radiation  trying  to  push  the  tran¬ 
sistor  into  conduction  for  V0  =  0  (i.e.,  decreases),  and 

the  radiation  causing  an  injection  of  negative  current  from 
the  diodes  into  the  base  of  the  transistor  for  V0  =  4  (resulting 
in  an  increase  of  V^) . 

The  overall  conclusion  from  the  gamma  testing  at  this 
level  is  that  the  XC-201  is  less  sensitive  to  radiation  than 
the  MC-201. 

c.  LINAC  Test  Results 

The  XC-201  and  MC-201  were  also  tested  on  the  LINAC  test 
facility  at  the  White  Sands  Missile  Range,  White  Sands,  New 
Mexico.  Both  of  the  units  irradiated  were  subjected  to  electron 
radiation  at  the  22-mev  electron  energy  level.  The  dose  measure¬ 
ments  for  the  LINAC  were  provided  by  WSMR  personnel  using  thermal 
luminescence  dosimetry.  There  was  considerable  variation  in  the 
indicated  dose  per  shot.  Average  figures  for  total  dose  per  shot 
were  about  120  rads  (l^O) .  The  dose  rate  may  be  calculated 
approximate ly  by  dividing  the  total  dose  by  the  radiation  pulse 
width.  For  the  electrons  on  the  LINAC  the  result  is  130  rads 

Q 

divided  by  0.5  Ms  or  approximately  2.6  x  10  rad/sec. 

In  figure  12,  the  external  test  circuit  and  conditions 
were  as  follows: 
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a) 

b) 

c) 

<*) 

e) 

f) 

The  data  were  recorded  by  a  Tektronix  type  555  oscilloscope. 
In  all  cases,  effectively,  the  microcircuits  only  were  irradiated. 
The  electron  beam  was  sufficiently  narrow  to  allow  the  positioning 
of  only  the  microcircuits  within  the  radiation. 

In  table  II  the  LINAC  results  are  summarized. 


Table  II 


LINAC  22 

Mev  Electron 

Results  on  the  MC 

-201  and  XC-201 

R^  Output 

Loading 

Resistance 

v8 

Input 

Voltage 

MC-201  Maximum 
Peak-to-Peak 
Transient  Volt¬ 
age  Output  (v^) 

XC-201  Maximum 
Peak-to-Peak 
Transient  Volt¬ 
age  Output  (Vj) 

Open 

Circuit 

0 

0.57  volt 

0.48  volt 

II 

1  volt 

no  shot 

0.20 

II 

1.5 

0.29 

no  shot 

II 

2 

no  shot 

0.16 

II 

5 

II 

0.10 

II 

0 

0.54 

0.52 

RL  =  1.6  K 

1 

no  shot 

0.24 

II 

2 

1.0 

0.14 

a  510-Cl  resistor  was  connected  from  terminals  6  to  5 
Vg  =  +8  volts 
Vc  =  -2.5  volts 

D 

Inputs  J,  9#  and  10  open 
Vg  variable 

The  effect  of  output  loading  is  shown  by  the  test 
for  equal  infinity  (i.e. 


,  open)  and  equal  1.6  K. 
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For  the  XC-201  the  transient  response  was  a  fast-rising 
negative-going  pulse  which  decayed  back  to  zero  with  a  time 
constant  of  approximately  3  x  10-^  sec.  The  peak-to-peak 
(i.e.,  here  the  same  as  the  peak)  maximum  transient  voltage 
output  (V-j)  is  shown  in  the  fourth  column  of  the  table.  Note 
that  the  output  response  rather  f  oi.^lusively  decreases  as  the 
input  Vg  increases.  This  is  as  to  be  expected  as  the  output 
transistor  is  cut  off  for  Vg  =  0,  but  is  going  toward  saturation 
for  Vg  increasing. 

For  the  MC-201  the  test  results  are  inconclusive.  Its 
transient  response  was  a  fast-rising  positive -going  pulse 
lasting  for  the  duration  of  the  radiation  pulse.  Then  the 
output  abruptly  switched  to  a  negative  pulse  which  also  decayed 
with  a  time  constant  of  approximately  3  x  10 sec.  Apparently 
two  simultaneous  phenomena  are  occurring,  one  pushing  the  output 
positive,  the  other  pushing  it  negative.  This  results  in  a 
cancellation  effect  which  for  low  values  of  input  voltage  makes 
the  peak-to-peak  response  of  the  MC-201  and  the  XC-201  approxi¬ 
mately  the  same.  Further  analysis  of  the  MC-201  will  require 
additional  shots  on  the  LINAC. 

d.  Summary  of  MC-201  and  XC-201  Tests 

In  general  the  XC-201  is  less  sensitive  than  the  MC-201  to 
transient  radiation  effects.  For  gamma  radiation  up  to  2  x  10^ 
roentgens/sec  these  results  are  conclusive.  For  electron  radia¬ 
tion  at  the  22-mev  level  the  difference  is  less  pronounced. 

3.  Fairchild  l.iL-905  Micrologic  Element  and  Mock-up 
a.  Background  Information 

The  Fairchild  micrologic  elements  are  a  set  of  compatible, 
integrated  logic  building  blocks.  The  elements  are  manufactured 
using  the  planar  epitaxial  process  by  which  all  the  necessary 
transistors  and  resistors  are  diffused  into  a  single  silicon 
wafer.  The  individual  RTL  gates  within  the  logic  elements  are 
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interconnected  by  metal  over  oxide.  These  units  are  character¬ 
ized  by  low  propagation  delays  for  use  in  high-speed  systems. 
Typical  propagation  delay  for  the  basic  RTL  circuit  is  12  nsec, 
and  pulse  widths  of  50  nsec  produce  reliable  triggers  for  the 
storage  elements. 

The  micrologic  gate  element,  mL-903#  is  a  three-input 
resistor-transistor-logic  (RTL)  circuit.  In  *  igure  13  are 
shown  the  circuit  and  logic  diagrams  for  the  dev*  **e.  For  a 
positive  input  on  1,  2,  or  3  the  associated  trax&aistor  goes 
into  saturation  which  reduces  the  potential  at  tne  output  (6) 
to  a  low  value  (i.e.,  the  "down"  sta  e).  For  no  positive 
inputs  on  1,  2,  or  3  the  gate  is  off  (i.e.,  the  "up"  state). 


Gate  Element 


(Top  View) 


Positive 
Logic : 


D  =  A+B+C 
=  A  B  C 


Negative 

Logic 


D  =  ABC 
=  A+B+C 


Figure  13 .  |-iL-903  Logic  and  Circuit  Diagrams 


b.  FXR  Tests  of  uL-903  and  Mock-up 


The  Kirtland  Air  Force  Base  flash  X-ray  unit  was  used  to 
irradiate  the  IJL-903,  The  tests  were  made  while  operating  the 
X-ray  tube  at  a  potential  of  400  kev.  The  incident  doses  were 

-7 

from  200  to  500  mi lliroentgens  with  a  pulse  duration  of  2  x  10 
sec.  The  dose  rates  were  of  the  order  of  10^  to  2.5  x  10^ 
roentgens/sec.  In  addition  to  irradiating  the  ijL-902  per  se,  a 
mock-up  of  the  device  was  constructed  from  discrete  diffused 
components  as  shown  in  figure  14.  The  components  of  the  mock-up 
are  similar  to  those  in  the  microcircuit,  each  diffused  on  a 
p-type  substrate.  The  mock-up  circuit  was  irradiated  in  whole 
and  in  part,  thereby  giving  3ome  insight  into  those  components 
which  showed  maximum  sensitivity  to  the  X-ray  pulse. 


Figure  14.  Mock-up  of  Fairchild  mL-903 
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c.  FXR  Test  Results 


Figure  15  shows  the  transient  radiation  response  of  the 
Fairchild  HL-905  circuit  and  its  mock-up  when  the  substrates 
(p-type)  of  the  circuits  were  grounded.  When  the  whole  circuits 
are  irradiated,  the  time  responses  of  both  circuits  are  almost 
identical.  The  sensitivity  of  the  mock-up  circuit  to  the  X  ray 
is  about  four  times  that  of  the  ML-905.  The  explanation  for 
this  lies  in  the  fact  that  the  active  volumes  of  the  transis¬ 
tors  in  the  mock-up  circuit  are  probably  greater. 

When  the  substrate  of  the  ijL-903  mock-up  is  tied  to  the 
source  of  the  collector  bias  voltage  (see  figure  15  h) ,  the 
diffusion  tail  shown  in  figure  15(a)  disappears.  The  collector- 
substrate  junction  is  in  a  forward  biased  state  ir.  this  case. 
Results  of  a  test  are  shown  in  figure  15(h) .  The  MET-1  transis¬ 
tor  was  in  the  cut-off  state  during  this  test  to  prevent  transis 
tor  burn-out,  because  this  arrangement  essentially  puts  V  on 
the  collector  of  the  transistor.  It  is  pointed  out  that  one 
does  not  normally  operate  a  circuit  with  the  substrate  tied  to 
V  .  This  was  only  a  diagnostic  procedure. 

When  the  substrate  of  the  nL-905  mock-up  is  tied  to  the 
collector,  the  collector  and  substrate  are  at  the  same  potential 
Under  operating  conditions  of  this  kind,  the  diffusion  tail  is 
about  1/5  of  its  value  when  tied  to  ground.  The  nET-1  was 
operated  in  the  cut-off  state  in  this  test.  The  test  results 
are  shown  in  figure  15(g). 

d.  Comments  on  FXR  Test  Results 

There  exists  a  negative  bias  at  the  substrate-collector 
junction  when  the  ML-905  is  operating  under  normal  conditions. 
When  this  condition  holds,  there  is  a  diffusion  tail  across 
the  barrier  with  a  time  constant  of  the  order  of  2.5  Msec. 

This  diffusion  tail  disappears  under  zero  bias  operation.  If, 
then,  zero  bias  operation  does  not  affect  circuit  performance, 
it  would  be  preferable,  from  a  point  of  view  of  radiation 
resistance,  to  operate  "nand"  elements  in  that  condition.  This 
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HL  905 


liL-903  Mock-^p 


(a)  D0B6  «  500  mr 
.002  V/div 
.5  usec/div 
inputs  Grounded 


(b)  Dose  w  200  mr 
0.1  V/dlv 
0.5  Msec/dlv 
inputs  Grounded 


(c )  Dose  m  500  mr 
.002  V/div 
0.1  Usec/div 
inputs  Grounded 


(d)  Dose  M  ?00  mr 
0.1  V/div 
0.1  Usec/cm 
Inputs  Grounded 


(e)  Dose  m  500  mr 
.002  V/div 


.002  V/div 
0.1  Usec/div 
Inputs  +1.5  Volts 


(f)  Dose  M  200  mr 
.005  V/div 
0.1  Msec/di  v 
inputs  +  1.5  Volts 


(9)  Dose  rate  «  200  mr 
.005  v/cm 
0.1  U sec/cm 
Inputs  Grounded 
Substrate  Tied  to  Collector 


(h)  Dose  rate  *,  200  mr 
.005  V/div 
0.1  wsec/div 
Input 8  Grounded 
Substrate  Tied  to  v 


Figure  15.  Transient  Response  of  the  uL-903 
"nand"  Gate  and  Its  Mock-up 


aq»an  points  out  the  importance  of  isolating  individual 
ocrHip^^onts  on  the  microcircuit  chip. 

The  identical  shape  of  the  responses  of  the  two  circuits 
j^Acates  that  this  mock-up  method  is  better  than  that  of 
using  conventional  components.  It  also  appears  to  be  a  method 
of  studying  more  complicated  microcircuits  in  a  radiation 
environment. 

e .  LIN AC  ML-903  and  Mock-up  Tests 

The  mL-903  and  its  mock-up  (see  figure  14)  were  irradiated 
at  the  WSMR  LINAC  facility.  The  4L-903  was  irradiated  by 
electrons  and  the  mock-up  was  irradiated  by  both  electrons  and 
gamma  rays.  All  of  the  tests  were  performed  at  the  22-mev 
electron  energy  level.  Average  figures  for  total  dose  are 
about  10  rads  (l^O)  for  gamma  and  about  130  rads  (l^O)  for 
electrons . 

In  table  III  the  results  of  the  HL-903  test  are  summarized. 
In  most  cases,  data  were  taken  with  no  load  (l^  =  °°),  fan-out 
of  one  (l^  =  500  O)  and  maximum  fan-out  of  five  (^  =  100  0). 

A  1N645  diode  was  additionally  inserted  in  series  with  the 
resistor  loads  to  better  simulate  the  ML-903  input  loading. 

Tests  during  the  irradiation  testing  at  the  flash  X  ray  indi¬ 
cated  that  these  resistor-diode  combinations  adequately 
represented  the  ml-903  inputs  as  loads. 


Table  III.  ML-903  Electron  Radiation- Induced  Response 


Output 

response 

Load 

<*L> 

For 

V,  =  0: 

For  Vj  =  1.5: 

00 

1.2 

volts 

1.2  volts 

«L  = 

500  a 

0.4 

0.52 

= 

100 

0.3 

0.28 

4l 


In  this  test  the  input  (v^)  was  etth#*  0  volt  or  1.5  volts. 
The  response  was  a  negative-going  pulse  lasting  only  slightly 
longer  than  the  radiation  pulse.  As  expected,  the  transient 
pulse  response  decreases  for  increased  loading. 

The  ML -903  mock-up  was  also  irradiated.  The  results  are 
summarized  in  table  IV.  Note  in  the  test  results  that  the 
collector  resistor  seems  to  be  the  most  sensitive  element. 

Two  gamma  ray  exposures  were  also  made.  The  output  peaks 
were  respectively  0.05  volt  and  0.04  volt  for  =  00  and 

-  500  O.  For  these  two  exposures  the  whole  mock-up  circuit 
was  irradiated. 

Table  IV.  LINAC  Electr  on  Test  of  ML-905  Mock-up 


Load 

Resistance 

0 

11 

K\ 

> 

=  1.5 

K  = 

GO 

0.4  volt  (collector 

0.15  volt  (input 

Jj 

resistor  only  irradiated) 

2  resistor  only) 

O.56  volt  (center  tran¬ 
sistor  only  irradiated) 

0.2  volt  (input  2 
resistor  only  irradiated) 

«L  “ 

500  O 

0.2  volt  center  tran¬ 
sistor  only  irradiated) 

No  data 

0.06  volt  (input  2 
resistor  only  irradiated) 

4.  TI  -  Series  51  Microcircuits 
a.  Introduction 

Some  of  the  Texas  Instruments  "Series  51"  microcircuits 
were  irradiated  with  the  AFWL  600-KV  flash  X-ray  machine. 
Tests  on  the  different  "Series  51"  circuits  yielded  similar 
results.  Therefore,  only  the  SN  512  "nand"  or  "nor"  gate 
was  studied  in  detail.  A  circuit  similar  to  the  SN  515  was 
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was  built  up  using  conventional  components.  The  transient 
radiation  effect  of  the  mock-up  circuit  compared  favorably 
with  the  SN  513  when  the  parasitic  transistors  were  taken 
into  account. 

b.  Test  Details 

The  Texas  Instrument  SN  512  and  SN  513  are  six-input 
"nor"  or  "nand"  logic  networks  made  from  diffused  silicon. 

The  schematic  of  the  SN  512  and  SN  513  is  shown  in  figure  16. 
Transistors  thru  are  input  transistors;  Qj  is  an  emitter- 
follower  output  transistor.  is  part  of  the  emitter-follower 

output  and  supplies  the  biasing  for  the  output  transistor. 

This  emitter-follower  output  circuit  is  useful  in  switching 
applications  where  a  low  collector  leakage  current  is  required 
and  a  negative  supply  voltage  is  not  available  for  reverse 
biasing  the  base  of  the  transistor.  A  six-volt  bias  was  used 
as  the  operating  voltage.  Normally  the  circuit  has  a  rise 
time  of  about  150  nsec  and  a  fall  time  of  1.4  l-isec.  The 
circuit  operates  in  the  current  mode,  which  means  it  is  always 
in  the  active  region  of  operation.  It  appears  that  in  the 
current  mode  of  operation  the  failure  mechanisms  caused  by 
radiation  may  be  different  for  the  "on"  state  and  "off"  state. 

When  the  circuit  is  operated  with  the  emitter-follower 
output  (SN-513)/  it  has  a  maximum  fan-out  capability  which 
will  drive  25  other  inputs. 

This  circuit  was  irradiated  in  the  flash  X-ray  beam  with 
each  exposure  delivering  about  500  mr  at  400  KV.  The  circuit 
was  operated  in  the  "off"  mode  by  connecting  all  inputs  to 
ground,  and  in  the  "on"  mode  by  applying  a  positive  1.5  volts 
to  the  input  of  the  base  circuit  for  each  transistor.  Both 
outputs  (8)  and  (9)  of  figure  16  were  monitored  by  using  a 
six-foot  length  of  RG-58  cable  terminated  in  a  series  combina¬ 
tion  of  a  51-tt  resistor  and  a  47-MMf  capacitor  at  the  moni¬ 
toring  oscilloscope.  The  capacitor  was  put  in  so  as  not  to 
interfere  with  the  DC  bias  level  of  the  circuit.  The  results 
are  shown  in  the  first  column  of  figure  17.  The  maximum 
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Figure  16.  Circuit  Diagram  of  SN  512  and 
SN  515  Micrologic  Circuit 


©  © 
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2-Resistor 

Microcircuit  Collector  Mock-up 


6-Resistor 
Diode  Mock-up 


» 

* 

• 

'  ^ ; 

i 

a  I 

i 

r 

V/div 

=  0.05  V/div  * 

0.005  V/div 

=  0.005 

All  inputs  are  grounded. 

Output  viewed  at  terminal  (8)  or  equivalent  of  figure  16. 


v/div  =  0.05 


V/div  =0.005 


V/div  =  0.002 


All  inputs  are  +1.5  volts. 

Output  viewed  at  terminal  (8)  or  equivalent  of  figure  16. 


All  inputs  are  grounded. 

Output  viewed  at  terminal  (9)  or  equivalent  of  figure  16. 


V/div  =0.005 


V/div  =0.002 


V/div  =0.002 


All  inputs  are  +1.5  volts. 

Output  viewed  at  terminal  (9)  or  equivalent  of  figure  16. 


Figure  17.  The  Response  of  SN  513  Microcircuit  and  Mock-up  to 
Flash  X-Ray  Pulses.  All  Time  Scales  are  0.1  usec/div.  The 
Dose  was  300  mr.  Vcc  =  6  Volts.  The  50  0  Cable  was  Terminated 
With  a  51  D  Resistor  in  Series  With  a  47  uuf  Capacitor 


45 


Figure  18.  Circuit  Representing  TI  SN  513  Made  of 
Conventional  Components 


Output  (8) 
Input  Grounded 


Input  -  +  1.5V 

All  sensitivities  ■  .005  V/div. 
All  time  scales  a  0.1  1 


t  i  i  >  >  i 


Output  (9) 
Input  Grounded 


*■  ^L>  i  *  <  I  i 


Output  (9) 
Input  =  +  1.5V 


Figure  19 .  Radiation  Pulses  of  SN  515  Mock-up 
Made  of  Conventional  Components 


fan-out  at  the  emitter-follower  output  (terminal  9)  is  25, 
and  at  the  collector  output  (terminal  8)  it  is  5.  Since  the 
input  resistance  of  similar  circuits  is  about  1000  0,  the 
range  of  driving  point  impedance  ranges  from  about  40  12  to 
1000  0. 

From  figure  17  it  can  be  seen  that  when  the  gate  transis¬ 
tors  are  in  the  "on"  state  the  transient  pulse  has  a  positive 
component  which  is  not  present  when  the  gate  transistors  are 
in  the  cut-off  state. 

An  SN  5135  circuit  was  built  using  conventional  transistors 
of  the  same  general  type,  carbon  resistors,  and  back  biased 
diodes  for  capacitors.  This  circuit  is  shown  in  figure  18. 

The  results  of  irradiating  the  circuit  are  shown  in  figure  19. 
The  radiation  effect  when  the  inputs  are  grounded  is  to  try 
to  turn  the  transistors  on.  There  was  a  5“miHivolt  peak  pulse 
at  output  (8) .  The  emitter-follower  output  under  the  same 
radiation  conditions  (500-mr  dose)  was  about  2  millivolts. 

When  the  inputs  were  biased  1.5  volts  positive,  a  similar 
effect  was  noticed,  indicating  that  the  transistor  tended  to 
be  driven  harder  into  saturation. 

At  this  point  it  was  obvious  that  there  was  something 
present  in  the  microcircuits  that  was  not  in  the  mock-up 
circuit  because  the  positive  peaks  did  not  occur  when  the 
gate  transistors  were  in  the  "on"  state.  The  next  step  was 
to  consider  the  parasitic  effects  in  the  microcircuits  and, 
if  possible,  to  introduce  them  into  the  mock-up  circuit. 

c.  Parasitic  Effects  in  SN  $12  and  SN  513 

In  figure  20  one  has  a  simple  transistor  circuit  which  is 
constructed  in  microcircuit  form  as  shown  in  figure  21.  Careful 
examination  of  figure  21  shows  that  for  both  the  transistor  and 
the  resistor  elements  there  exists  a  p-n-p  sequence  of  junctions 
starting  from  the  substrate  and  ending  at  the  base  of  the  n-p-n 
transistor. 


T 


i 


The  terminals  marked  N  and  S  correspond  to  the  isolation 

region  for  the  resistor  and  the  substrate.  N  is  usually  tied 

to  the  V  terminal  of  the  device,  and  S  is  tied  to  ground, 
cc 

In  this  case  the  parasitic  transistor  associated  with  the 
resistor  will  always  be  cut  off,  since  its  base  emitter  poten¬ 
tial  is  either  zero  or  reverse  biased.  The  parasitic  p-n-p 
transistor  associated  with  the  n-p-n  transistor  is  also  cut 
off,  provided  the  n-p-n  transistor  is  not  saturated.  The 
equivalent  circuit  is  shown  in  figure  22. 

Since  a  current  mode  switch  does  not  normally  saturate, 
it  is  not  necessary  to  resort  to  gold  doping  in  order  to  avoid 
parasitic  transistor  action.  Actually,  an  active  p-n-p  tran¬ 
sistor  action  will  exhibit  a  beneficial  effect  here,  in  that 
it  will  prevent  the  circuit  from  going  into  hard  saturation  if 
the  operation  of  the  circuit  is  such  that  it  is  overdriven. 

d .  Complete  Mock-up  of  the  SN  53-5 

To  simulate  the  parasitic  transistors  in  the  SN  5 15* 
some  2N995  p-n-p  transistors  were  put  in  the  circuit  as  shown 
in  figure  25*  The  collector  resistor  was  divided  into  two 
500-ohm  resistors  with  the  parasitic  transistor  emitter  connected 
between  the  two.  A  comparison  of  the  radiation  effects  between 
the  actual  microcircuit  and  the  more  complete  mock-up  is  shown 
in  figure  17.  Also  included  for  comparison  is  the  six-part 
collector  resistor,  back  biased  diode  model.  This  circuit  is 
shown  in  figure  24.  The  simple  two  collector  resistor  model 
produced  similar  radiation  effects  except  for  part  "e",  figure 
17.  This  positive  part  of  the  pulse  in  the  microcircuit  is 
probably  due  to  the  capacitance  associated  with  the  diffused 
collector  resistor.  The  six-resistor,  back  biased  diode  model 
had  too  much  capacitance  in  the  back  biased  junctions  to 
accurately  simulate  the  microcircuit  with  the  diffused  collector 
resistor.  A  model  somewhere  between  these  two  would  be  a  fair 
approximation . 
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Figure  25-  Mock-up  Circuit  of  SN  5 15  Using  Conventional 
Components  and  p-n-p  Transistors  to  Represent  Parasitic 
Transistors.  The  Collector  Resistor  is  Divided  Into 

Two  500  0  Resistors 
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Figure  24.  A  Six -Resist or  and  Diode  Conventional  Component 
Mock-up  Circuit  of  TI  Microcircuit  SN  512 
(Only  Two  Inputs  Shown) 
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e.  LINAC  Test  of  TI  SN  511  and  SN  5  HA 

The  Texas  Instruments  SN  $11  and  SN  5 HA  are  diffused  sil¬ 
icon  microelectronic  bistable  reset-set  flip-flop  counter  net¬ 
works  .  They  find  applications  in  digital  computer  systems , 
data  handling  systems ,  and  control  systems.  The  only  differ¬ 
ence  between  the  SN  511  and  SN  5HA  is  in  the  packaging. 

In  figure  25  is  shown  the  circuit  diagram  of  the  SN  5H- 
and  SN  5 11A.  These  circuits  were  irradiated  on  the  LINAC  at 

White  Sands  Missile  Range.  The  supply  voltage  used  was  V  = 

,  cc 

6  volts. 


preset  @  ®  CP 

Figure  -25.  TI  SN  5 11  and  SN  511A  Schematic 

The  tests  were  performed  with  the  output  10  in  the  down 
state  and  output  6  in  the  up  state.  The  results  for  both  gamma 
ray  irradiation  and  electron  irradiation  were  essentially  the 
same.  Output  10  was  a  positive  output  pulse  of  approximately 


52 


®  <§) 


7.5“Msec  duration.  After  7»5  Msec  the  output  10  returned  to  its 
original  down  state.  Output  6  was  a  negative  output  pulse  with 
an  exponential  decay  (tc  =  15  l^sec)  back  up  to  its  original  up 
state  level.  The  irradiation  tended  to  try  to  switch  the  flip- 
flop  but  it  restored  itself  to  its  original  state. 

The  tests  were  made  with  no  loading  on  the  circuits  other 
than  the  oscilloscopes.  In  table  V  the  peak  values  of  the  out¬ 
put  pulses  are  shown. 


Table  V 

SN  511  and  SN  5 HA  Radiation-Induced  Outputs 


Device 

Radiation 

Output  -  10  Peak 

Output  -  6  Peak 

SN  511 

Gamma 

+5-0  volts 

-2  .7  volts 

SN  5H 

Electrons 

+2  .0 

“2.7 

SN  511A 

Electrons 

+2.6 

«  -4 

5.  Signet ics  SU-315K 

a .  General  Information 

The  SU-315K  is  an  integrated  "nor"  circuit  fabricated 
within  a  monolithic  silicon  substrate  by  the  planar  technique. 
The  circuit  is  designed  for  use  in  high-speed  logic  applications 
and  is  designed  for  a  maximum  of  flexibility.  Its  input  and 
output  characteristics  are  compatible  with  those  of  other 
Signetic  SU-300  series  circuits .  Its  electrical  features  in¬ 
clude  high  noise  immunity,  high  speed  with  capacitive  loading, 
and  high  fan-out . 

The  SU-315K  is  a  dual  5-input  "nor"  gate.  Average  propa¬ 
gation  time  delays  run  from  20  to  40  nsec  depending  upon  the 
loading.  This  gate  has  a  maximum  fan-out  capability  of  12 
source  loads  plus  5  sink  loads. 
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Figure  26.  SU-315K  Dual  "Nor"  Gate  Schematic 


To  Amplifier 


+  4.5  v 


To  Amplifier 


(0  to+4v) 

(0  to  +  4v) 


Figure  27*  SU-315K  Tesc  Circuit 


b.  Test  Circuit 


Figure  26  shows  the  schematic  of  the  SU-315K  dual  "nor" 
gate.  There  are  two  identical  circuits  within  the  same  sub¬ 
strate.  The  half  of  the  circuit  associated  with  inputs  3#  4, 
and  5/  as  well  as  output  10 ,  was  actually  tested.  Inputs  3#  4, 
and  5  control  the  logic  output  state  cf  output  10.  is  the 
supply  voltage  for  the  transistor  collectors.  The  absolute 
maximum  voltage  rating  for  all  the  inputs  is  5  >5  volts. 

Figure  27  shows  the  actual  test  circuit.  Inputs  4  and  5 
were  grounded,  and  the  logic  state  of  output  10  was  controlled 
by  input  3«  The  logic  "1"  output  voltage  is  greater  than  3-3 
volts  for  V-j  ^  1.4  volts.  The  logic  "0"  output  voltage  is  a 
maximum  of  0.6  volt  for  V,  =  0.  Loading  impedances,  Z_  ,  were 

j  L 

placed  on  the  output  to  simulate  varying  fan-out  conditions. 

c .  FXR  Test  Results 

The  SU-315K  was  tested  on  the  AFWL  flash  X-ray  system. 

The  tests  were  performed  at  the  400-KV  level  with  a  0.2 -Msec 
pulse  width.  This  yields  an  approximate  dose  rate  of  2  x  10^ 
roentgens/sec . 

Table  VT  summarizes  the  SU-315K  FXR  test  results.  As 
shown  in  the  table,  shots  were  taken  for  varying  values  of  Vg, 
V-j ,  and  ZL . 

Shots  695#  690,  and  694,  in  order,  correspond  to  increas¬ 
ing  values  of  Z_ .  As  expected,  increasing  values  of  Z_  result 

li  ii 

in  an  increase  in  the  negative  peak  value  of  the  transient  out¬ 
put  and  an  increase  in  the  time  constant  of  the  decay. 

Shots  691#  697#  and  698  indicate  that  in  the  logic  "1"  state 
this  circuit  is  highly  radiation  resistant  to  gamma  dose  rates 
up  to  2  x  10^  roentgens/sec.  For  shot  692  the  SU-315K  was  set 
to  just  below  the  switching  level  (from  "0"  to  "1")  by  experi¬ 
mentally  adjusting  to  1.3  volts.  Irradiation  at  this  level 
produced  maximum  transient  output  for  a  fixed  load.  Shots 
690,  692,  and  691  indicate  that  the  transient  output  increases 
as  is  increased  up  to  the  point  of  switching.  Increasing 
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V-j  beyond  that  point  results  in  the  transient  output  being  re¬ 
duced  to  a  negligible  value.  Shots  69O  and  696  indicate  that 
the  transient  output  decreases  (within  certain  limits)  as  the 
supply  voltage  (Vg)  is  increased.  Shot  705  was  made  with  the 
device  behind  a  lead  brick.  As  this  shot  resulted  in  only  a 
small  output ,  it  can  be  concluded  that  the  stray  voltage  pick¬ 
up  contribution  to  the  output  pulses  is  correspondingly  small. 

Table  VI 

SU-315K  FXR  Test  Results 


Shot 

No. 

v6 

v3 

ZL 

Transient  Peak  of  Output  V^q 

690 

4 

0 

1.6  K 

ii 

75 

pf . 

0.02  volt 

691 

4 

4 

1.6  K 

11 

75 

pf. 

Negligible 

692 

4 

1.3 

1.6  K 

!! 

75 

pf . 

0.03  (V3  set  just  below 
switching  level) 

694 

4 

0 

1  meg. 

0.025  long  decay 

695 

4 

0 

500  n 

0.01 

696 

6 

0 

1.6  K 

I! 

75 

pf . 

0.015 

697 

6 

4 

1.6  K 

11 

75 

pf . 

Negligible 

698 

0 

4 

1.6  K 

1! 

75 

pf  • 

Negligible 

705 

0 

0 

1.6  K 

1! 

75 

pf . 

0.001  (device  behind  lead 
brick ) 

d.  LINAC  Test  Results 

The  SU-315K  was  also  tested  on  the  WSMR  LINAC  test  facility. 
The  unit  was  subjected  to  gamma  radiation  with  the  LINAC  at  the 
22-mev  electron  energy  level.  Average  figures  for  total  dose 
per  shot  were  about  10  rads  (l^O)  for  the  gamma  ray  mode.  The 
dose  rate  may  be  calculated  approximately  by  dividing  the  total 
dose  by  the  radiation  pulse  width.  For  the  gamma  LINAC  this 
result  is  10  rads  divided  by  0.5  nsec  or  approximately  2  x  10^ 
rads /sec . 
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Table  VII  gives  the  results  of  the  gamma  LINAC  testing. 

The  results  indicate  that  at  this  high  radiation  level  the  trans¬ 
ient  output  response  is  independent  of  the  logic  state  and  in¬ 
put  voltage  (v^).  Further  confirmation  of  this  conclusion  would 
be  interesting  if  additional  tests  are  scheduled  on  the  LINAC. 

Table  VII 

SU-215K  LINAC  Test  Results 


Shot 

No. 

v6 

v3 

ZL 

Transient  Peak  of  Output  V^q 

269 

4.5  volt 

0  volt 

Open 

0.05  volt 

270 

4.5 

2 

Open 

0.05 

271 

4.5 

3-5 

Open 

0.05 

e .  Summary  of  SU-515K  Testing 

The  SU-J15K  is  insensitive  to  gamma  radiation  while  in  the 
logic  "1"  state  up  to  at  least  2  x  roentgens/sec.  It  is 
sensitive  to  radiation  at  this  level  while  in  the  logic  "0" 
state . 

At  high  radiation  levels  the  SU-315K  becomes  sensitive  to 
gamma  radiation  irrespective  of  its  logic  state. 

6.  FXR  Test  of  Crysta Ionics  2N3084  FET 

The  2N3084  is  a  silicon  epitaxial  junction  N-channel  field 
effect  transistor  which  is  specifically  designed  for  use  in 
high-impedance  amplifiers  that  require  ultra-low  gate  leakage 
currents.  The  absolute  maximum  gate  leakage  is  0.1  nA.  It  is 
manufactured  by  the  epitaxial  junction  process  which  combines 
the  advantages  of  alloy,  epitaxial,  and  planar  techniques. 

The  device  is  bed  mounted,  oxide  passivated,  and  utilizes  a 
gold  bonding  process  which  eliminates  the  "purple  plague." 
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The  2N3084  was  tested  at  the  AFWL  flash  X-ray  facility. 
Both  the  gate  current  and  the  drain  current  radiation-induced 
transient  responses  were  recorded.  The  test  circuit  is  shown 
in  figure  28. 


IG  measurement 

^ - r —  f  VW 

!  1 

'  T _ 

Figure  28.  2N3084  Test  Circuit 

In  table  VIII  the  test  results  for  the  gate  current#  I_# 

VJ 

are  presented.  The  transient  response  was  a  negative  pulse 

r 

lasting  slightly  longer  than  the  radiation  pulse.  The  test  re 
suits  indicate  that  the  transient  peak  of  IG  increases  as  VQS 
decreases  and  that  the  peak  of  I_  increases  as  V__  decreases. 

V?  Do 

Table  VIII 

Peak  Value  of  Gate  Current  from  2N3084 


2N3084 
G 


D 

S 


I 


Ip  measurement 


VDS 

VGS 

Radiation 

I„  Peak 

V9 

10  volt 

1  volt 

yes 

0.013  ma 

10  volt 

2 

yes 

0.011 

10  volt 

5 

yes 

0.010 

10  volt 

4 

yes 

0.010 

2 

4 

no 

0.002 

2 

2 

yes 

0.013 

In  table  IX  the  test  results  for  the  drain  current#  Ip# 
are  presented.  The  transient  response  was  a  negative  pulse 
lasting  slightly  longer  than  the  radiation  pulse.  The  test 
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results  indicate  that  the  transient  peak  of  ID  increases  with 
increasing  VDg  at  the  VQg  =  4  volt  level,  and  ID  decreases  with 
increasing  VDg  at  VGg  =  0  volt  level. 

Table  IX 

Peak  Value  of  Drain  Current  (lD) 


VGS 

CO 

Radiation 

ID  Peak 

10  volt 

no 

0.0015  ma 

0  volt 

10 

yes 

0.0055 

15 

yes 

0.0035 

20 

yes 

0.0025 

4 

10 

yes 

0.002 

20 

yes 

0.004 

7 .  FXR  Test  of  Diffused  Resistor 

A  silicon,  phosphorous,  diffused  resistor  was  obtained  from 
the  microcircuit  laboratory  at  Rome  Air  Development  Center, 
Griffis  AFB,  New  York.  The  impurity  concentration  was  10  /cnr  , 
and  about  1.5  to  2  microns  deep.  The  resistance  was  4  ohms  per 
square  cm.  There  was  a  100-ohm  resistance  between  points  1  or 
2  and  3»  and  no  connection  to  the  substrate.  See  figure  29 • 

When  irradiating  this  resistor  with  a  dose  rate  of  5  x  10^ 
r/sec,  or  a  total  dose  of  about  1  r  in  air,  a  peak  transient 
current  of  l6  microamperes  was  obtained.  The  current  pulse  fol¬ 
lowed  the  X-ray  pulse  except  for  a  tail.  This  tail  was  init¬ 
ially  2  microamperes  and  decayed  out  in  about  5  nsec.  The  same 
pulse  was  observed  when  the  resistor  was  irradiated  in  oil  to 
eliminate  any  possible  air  effects.  The  oscilloscope  trace  is 
shown  in  figure  50.  The  diffused  resistors  appear  to  have  a 
large  transient  current  similar  to  that  observed  in  microcir¬ 
cuits;  the  current  pulse  is  probably  caused  by  the  isolation 
junction. 
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Figure  29*  Silicon  Diffused  Resistor 


Figure  JO.  Transient  Radiation  Response 
of  a  Diffused  Silicon  Resistor 
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3.  Thin-Film  Nichrome  Resistors 

Several  thin-film  nichrome  resistors  ranging  in  resistance 
value  from  16  ohms  to  11.4  K  were  irradiated  in  the  600-KV  flash 
X  ray.  No  effect  was  observed  on  the  smaller  resistance  values 
in  a  radiation  field  up  to  5  x  10^  r/sec.  The  11.4-K  resistor 
has  a  radiation  response  in  air  that  is  similar  to  a  carbon 
resistor.  The  response  is  due  to  air  ionization  effects.  When 
the  11.4-K  nichrome  resistor  was  irradiated  in  oil,  most  of  the 
effect  observed  in  air  disappeared.  There  was  a  small  pulsed 
response  of  about  2 .5  microamp  peak  and  1  microsecond  in  dura¬ 
tion.  This  is  (most  likely)  secondary  emission  of  electrons  and 
the  resulting  replacement  current. 

9.  The  WSMR  Gamma  LINAC  Experimental  Results  on  Hall  Devices 

a .  General  Results 

The  theory  of  the  operation  of  Hall  devices  in  a  radiation 
^.ivironment  was  developed  and  reported  in  AFWL  Report  No.  TDR- 
64-58.  However,  the  experimental  data  was  not  complete  until 
the  present  WSMR  gamma  LINAC  was  made  available. 

The  theory  predicted  that  there  should  be  no  radiation  ef¬ 
fects  for  the  constant  voltage  mode  of  operation  and  that  the 
Hall  voltage  should  vary  as  a  function  l/n  for  the  constant  cur¬ 
rent  bias  condition,  where  n  is  the  number  of  free  carriers  per 
unit  volume  of  the  Hall  material. 

The  thin-film  Hall  devices  were  found  to  be  highly  radia¬ 
tion  resistant  as  was  predicted  in  the  theory.  There  were  no 
radiation  effects  observed  for  dose  rates  up  to  10^  r/sec 
when  the  devices  were  operated  in  the  constant  voltage  mode. 

This  is  the  operating  mode  that  should  be  used  for  any  Hall  de¬ 
vice  where  gamma  radiation  environments  or  extreme  temperature 
changes  might  be  encountered. 

When  operating  Hall  devices  in  the  constant  current  mode, 
the  observed  effects  followed  closely  those  predicted  in  the 
above  report.  The  percentage  change  in  Hall  voltage  for  the 
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constant  current  case  is  given  in  table  X.  The  results  from  the 
LINAC  experiments  are  shown  in  figures  31  and  32  for  the  HS-51 
and  figures  26,  27#  and  33  for  the  Beckman  device. 


Table  X 

Hall  Voltage  Change  (Constant  Current)  During  Irradiation 


Dose  Rate 
• 

y 

Roentgen/Sec 

Avh* 

Device 

HS-$1 

Device  335 

Device  331 

exp. 

theory 

exp. 

theory 

exp. 

theory 

108 

0 

9.1  x  10“8 

0 

- 

0 

- 

io9 

0 

.01 

0 

- 

0 

- 

1010 

5 .9  to  9 

9 

0.75 

0.8 

0 

0.187 

1011 

- 

50 

- 

7.5 

- 

3.44 

1012 

- 

90 

44.8 

- 

15-76 

1013 

- 

99 

- 

89 

- 

65.28 

The  Hall  voltage  change  in  a  pulsed  radiation  environment 
for  the  constant  current  mode  of  operation  is  given  in  fig¬ 
ure  3 Id.  The  percentage  change  in  Hall  voltage  was  5.5  percent 
at  a  dose  rate  of  8.2  x  109  r/sec.  When  the  magnetic  field  was 
removed,  no  change  was  observed  as  can  be  seen  in  figure  31b. 

And  when  the  magnetic  field  is  reversed,  tne  AV^  (change  in 
magnitude  of  Hall  voltage)  is  in  the  opposite  direction  as  pre¬ 
dicted.  In  figure  31c  the  percentage  change  was  9  percent,  for 
a  dose  rate  of  9*2  x  10  .  If  the  bias  current  is  reversed,  the 
direction  of  the  Hall  voltage  is  also  reversed  and  a  similar  re¬ 
sult  was  obtained  (figure  32a).  The  theory  predicted  no  change 
in  Hall  voltage  for  a  constant  voltage  operating  mode  and  this 
is  what  was  observed  as  can  be  seen  in  figure  32b.  Referring 
to  the  constant  bias  current  case,  theory  demands  that  a  change 
of  input  impedance  equal  the  change  in  Hall  voltage.  This  was 
observed  and  the  results  are  shown  in  figure  32c.  The  radiation 
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a. 

B  *  4  kilogauss  N-S 
IB=  300  nia  (+  -) 

Y  =  8.2  x  109  r/sec 
t  -  2  Msec/cm 
^fivh  =  5 


b. 

B  =  0 
Ifi=  300  ma 

y  *  8.1  x  10?  r/sec 
t  =  2  ^  sec/cm 

»Vh  =0 


Radiation 


c . 


4  kilogauss 
300  ma 

9.2  x  109  r/sec 

2  |i  sec/cm 

s  9* 


Figure  31.  Radiation  Effect  on  the  Ohio  Semiconductor 

HS-51  Hall  Device 


a . 

B  =  4  kilogauss  S-N 
IB=  300  ma  (-  +) 

Y  =  9  x  109  r/sec 
t  =  2  Msec/cm 
#Wh  =  9 % 


b. 

B  ■  4  kilogauss  S-N 
VB=  0.3  Volts 

Y  *  9  x  109  r/sec 
t  *  2  |i  sec/cm 

4Vh  - 0 


c . 

B  =  4  kilogauss  S-N 
Y  -  9  x  109  r/sec 
6zin  =  IQ* 


Figure  32  .  Radiation  Effects  on  the  Ohio  Semiconductor 

HS-51  Hall  Device 


B  =  4  kilogauss 
I0=  20  ma 

Air  Gap  =0.7  cm 
#fiVh  =  1.37# 

Y  -  10 10  r/sec 


(10  MV)  \L 
cm 


Radiation 

Pulse 


b.  Constant  Voltage  Bias 


Radiation 

Pulse 


4Vh  =° 

Y  ■  101®  r/sec 


(10  MV)  V. 
cm 


c . 
6IB 


Input  Current 
-  150  ua 
or 

0.75# 


Radiation 

Pulse 


(500  iia)  IB 

cm 


Figure  33.  Radiation 


pulse  as  monitored  by  a  photodiode  and  plastic  scintillator  is 
also  shown  on  each  picture. 

It  was  found  when  irradiating  the  thin-film  devices  that  a 
much  larger  effect  is  observed  when  a  small  air  gap  is  used. 

This  indicates  that  secondary  electrons  knocked  out  of  the  mag¬ 
net  have  an  effect  on  the  Hall  device.  If  the  oscillographs  for 
the  Beckman  335  device  are  observed,  a  28-MV  change  is  recorded 
at  the  Hall  terminals.  In  this  case,  the  device  was  fitted 
closely  in  the  air  gap  of  the  magnet.  When  the  device  was 
mounted  in  the  center  of  a  three-centimeter  air  gap,  the  AV^ 
decreased  to  1.5-MV  change  or  0.75-P®rcent  change,  which  agreed 
with  the  values  predicted  from  theory  (  Table  X  ).  The  thin- 
film  devices  also  have  a  larger  input  impedance  and  the  current 
did  not  remain  constant  during  the  pulse.  The  impedance  of  the 
current  source  was  not  quite  large  enough  to  hold  the  current 
constant.  As  can  be  seen  from  figure  33c,  there  was  a  change 
in  input  current  of  about  150  microamperes  during  the  radiation 
pulse.  Figure  33a  shows  the  transient  effect  when  an  air  gap 
of  0.7  centimeter  was  used.  The  change  in  Hall  voltage  is  1.37 
percent  in  this  case  and  is  higher  than  the  predicted  value  of 
0.8  percent  indicating  that  secondary  electrons  are  present, 
but  not  as  much  as  in  the  smaller  air  gap.  Figure  33b  indicates 
no  change  in  Hall  voltage  for  the  constant  voltage  operating 
mode  which  is  the  same  as  observed  for  the  HS-51. 

Similar  results  were  obtained  using  the  Beckman  331  device. 
This  is  a  device  which  is  five  microns  thick,  and  the  predicted 
AV^  at  1010  r/sec  is  less  than  0.19  percent.  There  was  no 
change  observed  for  this  device. 

b.  Conclusions 

A  theory  was  developed  that  predicts  the  transient  radia¬ 
tion  effects  on  Hall  devices.  These  effects  depend  on  the  type 
of  biasing  used.  Experimental  results  agree  with  the  theory 
for  both  constant  voltage  and  constant  current  modes  of  opera¬ 
tion.  The  theory  for  the  constant  voltage  case  predicted  no 
change  in  Hall  voltage  and  none  was  observed.  The  highest  dose 
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rate  available  was  10^  r/sec.  Therefore,,  whenever  high- 
intensity  qamma  radiation  is  a  factor  in  the  environment ,  the 
constant  voltage  mode  of  operation  should  be  used.  This  would 
also  be  true  for  any  other  factor  in  the  environment  which  would 
cause  a  change  in  carrier  density,  such  as  large  temperature 
variations  . 

The  Hall  voltage  was  found  to  vary  as  l/n  for  the  constant 
current  mode  of  operation.  At  a  dose  rate  of  lO1^  r/sec  a  max¬ 
imum  change  in  Hall  voltage  of  9  percent  was  predicted  for  the 
Ohio  semiconductor  HS-51  device.  A  change  of  5*5  percent  to  9 
percent  was  observed. 

The  high-impedance  thin-film  devices  are  better  suited  for 
constant  voltage  drive  than  the  bulk-type  devices.  They  are 
also  more  radiation  resistant.  At  a  dose  rate  of  lO1^  r/sec 
they  had  a  percentage  change  in  Hall  voltage  of  at  least  one 
order  of  magnitude  less  than  the  bulk-type  devices  in  the  con¬ 
stant  current  mode  of  operation. 

10.  Summary  of  Conclusions 

Several  conclusions  may  be  drawn  from  this  research: 

(1)  In  general,  the  EPIC  process  XC-201  is  less  sensitive 
than  the  MC-201  to  transient  radiation  effects.  For  gamma  radi¬ 
ation  up  to  2  x  106  roentgens/sec  these  results  are  conclusive. 
For  electron  radiation  at  the  22-mev  level  the  difference  is 
less  pronounced. 

(2)  Testing  of  the  pL-903  indicated  that  the  substrate- 
collector  junction  introduced  a  diffusion  tail  in  the  transient 
radiation  response  which  can  be  reduced  by  elimination  of  the 
isolation  junction.  The  successful  mocking  up  of  the  mL-903 
showed  this  technique  to  be  a  useful  one  for  research. 

(3)  Mock-up  circuits  for  the  SN  513  required  the  inclu¬ 
sion  of  the  parasitic  elements  for  accurate  simulation. 
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(4)  The  SN  $11  and  SN  511A  flip-flops  were  irradiated. 

The  irradiation  tended  to  try  to  switch  the  flip-flop,  but  it 
restored  itself  to  its  original  state. 

(5)  The  SU-315K  is  insensitive  to  gamma  radiation  while 
in  the  logic  "1"  state--at  least  up  to  2  x  10^  roentgens/sec. 

It  is  sensitive  to  radiation  at  this  level  while  in  the  logic 
"0"  state.  At  high  radiation  levels  the  SU-315K  becomes  sensi¬ 
tive  to  gamma  radiation  irrespective  of  its  logic  state. 

(6)  The  Crystalonics  2N3084  field  effect  transistor  was 
tested  on  the  600-KV  flash  X  ray.  The  test  results  indicate 
that  the  transient  peak  of  I.  increases  as  V__  decreases  and 
that  the  peak  of  I-  increases  as  V_0  decreases.  The  test  re- 
suits  also  indicate  that  the  transient  peak  of  ID  decreases  with 
increasing  VDg  around  the  Vgg  =  0  volt  level,  and  that  ID  in¬ 
creases  with  increasing  V__  at  the  V_„  =  4  volt  level. 

Do  Go 

(7)  Diffused  resistors  have  a  large  transient  current  sim¬ 
ilar  to  that  observed  in  many  microcircuits.  The  current  pulse 
is  probably  caused  by  the  isolation  junction. 

(8)  The  transient  radiation  response  of  thin-film  nichrome 
resistors  is  primarily  due  to  air  ionization  effects. 

(9)  Thin-film  Hall  devices  were  found  to  be  highly  radia¬ 
tion  resistant  as  was  predicted  by  the  theory.  There  were  no 
radiation  effects  observed  for  dose  rates  up  to  10^  r/sec  when 
the  devices  were  operated  in  the  constant  voltage  mode.  This 
is  the  operating  mode  that  should  be  used  for  any  Hall  devices 
where  gamma  radiation  environments  or  extreme  temperature  changes 
might  be  encountered. 
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SECTION  III 


PRELIMINARY  DISCUSSION  OF  ENGINEERING  MODELS 
FOR  THIN-LAYERED  DEVICES 

1 .  Introduction  to  Fundamentals  of  Diffusion  Theory 

The  purpose  of  this  discussion  is  to  establish  very  gen¬ 
erally  the  basis  of  the  diffusion  theory  used  in  semiconductor 
device  analysis.  Diffusion  theory  is  obtained  from  the  more 
rigorous  transport  theory  by  making  simplifying  approximations 
to  the  Boltzmann  transport  equation.  This  procedure  leads  to 
the  diffusion,  or  continuity,  equation  for  holes  and  electrons 
in  semiconductors.  The  assumptions  made  in  the  development  of 
the  theory  are  then  shown  to  be  false  in  their  application  to 
thin- layered  devices. 

2 .  The  Distribution  Function 

The  various  electronic  transport  properties  of  a  conductor, 
or  semiconductor,  can  be  found  if  the  behavior  of  the  electronic 
distribution  function  can  be  determined  under  the  influence  of 
externally  applied  fields.  The  electronic  distribution  function, 
f(r,  k;,  t),  gives  the  probability  at  time  t  of  the  occupation  of 
the  energy  state  corresponding  to  the  wave  vector  k  at  a  point 
in  the  crystal  given  by  the  position  vector  r.  More  rigorously, 
since  f  is  a  continuous  probability  density  function,  one  should 
say  that  the  quantity 

f(r,  k,  t)dr  dk  (23) 

is  the  probability  at  time  t  of  particles  occupying  the  volume 
element  dr  which  is  located  by  the  position  vector  r  and  in  the 
element  of  wave  vector  space  dk  about  the  wave  vector  k.  Using 
either  definition,  it  should  be  clear  that  the  function  f  gives 
one  the  information  necessary  to  determine  the  manner  in  which 
particles  change  energy  and  position  within  the  crystal. 


3 .  The  Boltzmann  Transport  Equation 


The  basic  problem  in  the  study  of  conductors  and  semicon¬ 
ductors  is  the  determination  of  the  function  f  under  the  influ¬ 
ence  of  externally  applied  fields  (electric,  magnetic,  and 
thermal).  This  problem  is  usually  attacked  through  the  Boltz¬ 
mann  transport  equation  which  described  the  manner  in  which  the 
particles  behave  in  the  material.*  Fundamentally,  the  transport 
equation  postulates  that  the  total  time  rate  of  change  of 
f(r,  k,  t)  is  caused  by  particle  collisions  with  the  fixed  atoms 
of  the  lattice.  Mathematically  this  is  stated  as 


df  _  df 
dt  "5tT 


coll. 


(24) 


Before  considering  the  collision  terms,  one  should  look  at  the 
total  derivative  of  f  given  by 


df  _  8f  dx 

at  "  ■sr  at 


.3fdy.3fdz.3f  dkx  .  3f  dky 
+  i73t+'5?at+'5!rd^  +  ^rdt 

x  y 


3f  dkz  .  3f 

w  at-  +'5t 

z 


(25) 


where  r  and  k  have  been  replaced  by  their  components  in  rec¬ 
tangular  coordinates.  Interpreting  dx/dt ,  dy/dt ,  and  dz/dt  as 
components  vx,  v^,  and  vz ,  respectively,  of  the  particle  veloc¬ 
ity  v ,  and  dkx/dt ,  dky/dt ,  and  dkz/dt  as  components  of  a  vector 
dk/dt  one  can  write  equation  (25)  as 


df  _ 
dt  " 


v 


v-f  + 
r 


0.  ■ 

dt 


vi?f 


(26) 


*For  more  detail  see  the  books  by  Wilson  (reference  20), 
Sietz  (reference  21),  and  Smith  (reference  22),  and  the  article 
by  Blatt  (reference  23 )• 
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where 


v-f 


r 


Bf  ,  Bf  -•  ,  Bf  -• 

Bx  ax  By  ay  Bz  az 


„  -  _  Bf  -*  ,  Bf  -  .  Bf  - 

k  Bk  x  Bk  y  Bk  z 

a  y  4 

with  a  ,  a ,  3nd  3  unit  vectors  in  the  respective  directions. 

a  y 

Since  electrons  in  3  crystal  behave  as  free  particles  with 
energy 


where  p  is  the  momentum  of  the  particle  and  m  the  effective 
mass,  one  can  transform  (26).  The  momentum  is  given  by  p  =  mv ? 
thus  dp/dt  =  d/dt  (mv)  and  if  m  is  a  constant  in  time, 

if  =  raiF  =  ma  =  r  (28) 

where  F  is  the  force  on  the  particle  and  a  is  its  acceleration. 
Bragg's  relation  relates  the  momentum  to  the  wave  vector  as 
p  =  h  k,  where  h  =  Planck's  constant/27r.  Thus, 

||  =  h  ^  =  F  implying  that  =  £  (29) 

and  using  (20)  in  (26)  one  finally  arrives  at 

||  =  v  •  v-f  +  £  F  •  Vjjf  +  ||  (50) 


The  collision  terms  to  which  (50)  must  be  equal  should  now  be 
considered.  The  collision  term  is  described  by  a  scattering 
process  which  takes  an  electron  from  a  given  state  it  to  a  state 
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ic'  and  vice  versa.  Let  w(ic,  k '  ,  r,  t)*  be  the  probability  that 
a  particle  in  about  r  at  time  t  will  change  its  wave  vector 
from  £  to  1c'  due  to  some  type  of  collision.  Then  to  include  the 
exclusion  principle  one  may  multiply  w(k,  k' )  by  the  probability 
that  the  initial  state  £  is  occupied,  f(£),  and  the  probability 
that  the  final  state  is  empty  (l  -  f(k')).  Then  the  probability 
of  a  transition  from  state  k  to  state  £'  becomes 

w(i?,  ic •  )f(i?)(i  -  f(i?' )).  (51) 

Equation  (51)  represents  a  loss,  by  collision,  of  particles  from 
an  original  state.  The  gain,  from  collisions,  to  the  original 
state  is  given  similarly  by 

w(ic\  ic)f (k'  )(1  -  f(ic)).  (52) 

The  total  change  in  f,  due  to  collisions,  at  r  and  k  is  given 
by  integrating  the  difference  of  (51)  and  (52)  over  all  k' 

The  collision  term  then  becomes 

(|f)  =  \[-w(S.  i?')f(Od  -  f(i?'))  + 

coll  J 

w(£,  ic')f(k')(i  -  f(£))]d£\  (55) 

Equating  (50)  to  (55)  gives  the  fundamental  equation  for  deter¬ 
mining  the  properties  of  conductors  and  semiconductors: 

|i  +  7  •  7?f  +  £  F  •  V£f  =  $[-w(S,  it'  )f(ic)(i  -  f(ic'))  + 

W(k,  lc •  )f (lc ■ ) ( 1  -  f(i?))]dk'.  (54) 


*Thc  symbol  W(k,  £' ,  r,  t)  will  be  labeled  w(k‘l  ,  k)  in  the 
interests  of  brevity. 
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Some  discussion  of  equation  (34)  is  in  order.  It  is  the 
type  of  classical  integro-dif ferent ial  equation  first  developed 
in  the  study  of  gas  dynamics.  Solutions  obtained  to  date  for 
this  equation  are  at  best  approximate  and  usually  done  only  in 
the  steady  state,  i.e.,  when  df/dt  =  0.  It  should  be  realized 
that  (34)  itself  is  only  an  approximate  relation  for  the  phe¬ 
nomena  that  occur  within  the  crystal.  It  is  the  best  known  and 
probably  the  most  widely  used  relation  in  the  study  of  metals 
and  semiconductors.  It  is  not  the  purpose  here  to  develop  this 
equation  but  merely  to  recognize  it  as  the  starting  point  for  the 
study  of  various  phenomena  in  metals  and  semiconductors.  The 
first  term  on  the  left  is  the  time  variation  of  the  density  func¬ 
tion  f.  The  second  term  introduces  the  effects  of  particle  den¬ 
sity  and  thermal  gradients  in  the  material.  The  third  term  on 
the  left  brings  in  the  effects  of  applied  electrical  and  magnetic 
fields  through  the  force  F,  which  in  the  case  of  charged  parti¬ 
cles  is  given  by 

F  =  q(E  +  V  x  B)  (MKS) 

where 

q  =  particle  charge 

E  =  applied  magnetic  field 

v  =  particle  velocity 

and 

B  =  applied  magnetic  field. 

Thus  the  left-hand  side  introduces  the  effects  of  the  externally 
applied  fields  and  of  particle  density  gradients.  The  right- 
hand  side  of  equation  (34)  is  at  best  a  guess  at  what  goes  on 
inside  the  material  on  a  microscopic  scale,  integrated  to  give 
the  macroscopic  effect . 

The  solution  to  equation  (34),  when  undertaken,  is  done  by 
approximate  means.  The  assumption  that  the  density  function  for 
the  disturbed  system,  f(k,  r,  t),  is  different  from  some  equi¬ 
librium  distribution  fQ(r,  £,  t)  by  only  a  small  amount  is 
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usually  made.  Such  an  assumption  allows  one  to  expand  f  in  a 
power  series  about  f  and  keep  only  the  first-order  expansion 
term.  Surprisingly  enough,  this  has  been  a  fairly  satisfactory 
approach  in  providing  agreement  with  experiments.  The  results 
of  this  approach  will  not  be  studied,  but  a  simplifying  assump¬ 
tion  will  be  made  which  shall  be  discussed  in  detail,  that  of 
the  existence  of  a  characteristic  relaxation  time  for  the 
material . 

a .  Relaxation  Time  and  the  Collision  Integral 

By  far  the  most  often  made  simplifying  assumption  in  dis¬ 
cussing  equation  (34)  is  to  represent  the  collision  integral  by 
a  process  of  relaxation  given  by 


(M) 


coll 


f  (ic,  r  ,  t )  -  fQ(£,  r ,  t ) 
T~ . . 


(56) 


Equation  (36)  implies  that  for  any  suddenly  applied  or  removed 

fields,  the  density  function  f  relaxes  to  an  equilibrium  density 

function  f  with  a  characteristic  time  t.  The  existence  of  a 
o 

relaxati<  a  time  immensely  simplifies  the  transport  problems. 

The  validity  and  consequences  of  the  assumption  of  a  relaxation 
time  shall  now  be  discussed  at  some  length. 

Whether  or  not  a  relaxation  time  can  be  defined  for  the 
collision  integral  depends  basically  on  two  things,  the  first 
being  the  assumed  nature  of  the  colxision,  and  the  second  being 
the  temperature  at  which  the  process  takes  place.  There  are 
generally  two  scattering  processes  encountered  by  the  particles, 
one  being  due  to  lattice  motion,  and  the  other  being  due  to  im¬ 
purities  and  lattice  defects.  The  lattice  scattering  can  be 
further  divided  into  particle  interactions  with  acoustical  or 
optical  modes  of  the  lattice  vibration.  In  the  interaction 
with  the  acoustic  mode  of  lattice  vibration,  the  energy  change 
of  the  charged  carrier  is  small  compared  to  its  average  energy. 
This  is  essentially  due  to  the  low-frequency  nature  of  the 
acoustic  phonon.  In  this  case  the  phonon  energy  ha)  «  kT  for 
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any  temperature  greater  than  the  Debye  temperature  of  the  mate¬ 
rial.*  Consequently,  during  a  phonon-particle  collision,  the 
particle  can  only  gain  small  amounts  of  energy.  Therefore, the 
density  function  f  will  not  differ  appreciably  from  the  equilib¬ 
rium  function  f  ,  and  the  collision  then  can  be  approximated  by 
the  relaxation  time  assumption. 

The  optical  mode  scattering  of  particles,  because  of  the 
higher  frequency  of  the  optical  phonons,  allows  the  relaxation 
time  approximation  only  at  higher  temperatures.  At  very  low 
temperatures  one  can  have  kT  <  he*,  and  a  phonon -particle  colli¬ 
sion  would  result  in  a  large  particle  energy  change.  At  temper¬ 
atures  where  ha)  «  kT,  the  same  condition  holds  as  for  acoustic 
phonons,  and  a  relaxation  time  can  be  defined.  It  has  also  been 
shown  that  for  nonpolar  solids,  a  relaxation  time  always  exists 
(reference  24).  Therefore, one  concludes  that  for  scattering  by 
lattice  vibrations  the  relaxation  time  can  be  defined  for  temp¬ 
eratures  in  excess  of  the  Debye  temperature. 

In  considering  scattering  by  impurities  and  lattice  defects, 
one  assumes  the  collisions  to  be  perfectly  elastic.  In  such 
cases  f(k')  =  f(k)  and  a  relaxation  time  can  always  be  defined. 

It  is  concluded  that  one  can,  in  general,  define  relaxation 
times  for  the  scattering  processes  if  the  energy  of  the  particles 
is  changed  only  a  small  amount  by  the  interactions.  The  first 
tv  due  to  lattice  vibration  is  temperature  dependent;  the  sec¬ 
ond  tt  due  to  impurities  and  defects  is  independent  of  tempera- 
ture.  These  are  combined  to  define  a  total  relaxation  time  t  by 


Now  one  may  consider  equation  (33)  and  see  under  what  con¬ 
ditions  a  relaxation  time  can  be  a  valid  representation  for  the 


*k  here  is  the  Boltzmann  constant  and  should  not  be  con¬ 
fused  with  the  wave  vector  k. 
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The  relaxation  time  comes  from  the  solution  of  this  linear  in¬ 
homogeneous  integral  equation  in  terms  of  f(k)  -  f  (k*)  with 
df/dt  =  0.  The  relaxation  times  which  can  be  obtained  depend  on 
the  form  of  the  left-hand  side  of  equation  (42). 

b.  Relaxation  Times  from  the  General  Solution  of 
Boltzmann's  Equation 

The  Boltzmann  equation  (42)  has  an  exact  analytic  solution 
under  certain  simplifying  circumstances.  The  first  of  these  is 
that  the  energy  surfaces  in  k-space  are  spherical  so  that  v  and 
£  point  in  the  same  direction.  This  allows  one  to  write  the 
left-hand  side  of  (42)  as  a  scalar  product  of  v  (with  a  vector 
which  is  constant  on  a  spherical  surface  in  k-space  and  is 
therefore  proportional  to  cos  ,  where  ip  is  the  angle  between  k 
and  the  constant  vector).  The  second  condition  is  that  the 
scattering  probability  depends  only  on  the  angle  between  k  and 
k'  and  not  on  their  absolute  orientations  in  the  crystal.  Then 
one  can  define  a  differential  transition  probability  o)(jc,  0)dfl, 
which  is  the  probability  that  a  particle  of  momentum  hk  is 
scattered  through  the  angle  0  into  the  solid  angle  dfl. 

Under  these  assumptions  it  can  be  shown  that  equation  (42) 
has  the  exact  solution  given  by 


•  Lll—  ?i<fo ^ l 

\(l  -  cos  9)co(ic,  0)dO 


(4?) 


Now  the  integral  in  the  denominator  of  (43)  is,  dimensionally, 
a  reciprocal  transition  rate  for  electrons  of  momentum  tik. 
Recognizing  this ,  one  writes 


— - —  =  ^(l  -  cos  0 )a> (ic ,  0)dO  (44) 

t(*) 

From  (44),  for  experimentally  determined  o>(k,  0),  a  valid  relax¬ 
ation  time  can  be  determined  which  provides  one  with  a  di~tribu- 
tion  function  f(£)  from  which  one  can  determine  the  various 
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electronic  transport  properties  of  the  material.  On  the  other 
hand,  from  measurements  of  thermal  and  electrical  transport 
characteristics  a  relaxation  time  can  be  determined. 

Examples  of  such  relaxation  times  determined  from  measure¬ 
ments  of  electrical  and  thermal  conductivities  are: 


electrical 


(45) 


thermal 


(46) 


In  these  last  relations  the  notation  is: 

m  =  effective  mass 
n  =  density  of  electrons 
e  =  electronic  charge 
k  =  Boltzmann  constant 
o  =  electrical  conductivity 
K  =  thermal  conductivity 
T  =  absolute  temperature 

In  general  Te  and  are  not  equal  except  in  cases  where  the 
Wiedemann-Franz  Law  holds.  Thus  the  fact  that  the  Lorenz  num¬ 
ber  L,  defined  by 


=  (*7) 

is  a  constant  independent  of  the  properties  of  the  material  is 
a  consequencesof  the  existence  of  a  unique  time  of  relaxation. 
The  fact  that  at  temperatures  below  the  Debye  temperature  L  is 
no  longer  a  constant  is  due  to  the  relatively  large  energy  ex¬ 
changed  in  electron-phonon  interactions.* 


*For  a  more  detailed  discussion  see  (reference  2J>,  pp. 
222-224)  or  (reference  20,  pp.  201-202). 
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It  is  clear  from  (44)  that  the  type  of  scattering  process 
will  determine  the  relaxation  time.  Equation  (43)  need  only  be 
satisfied  on  each  constant  energy  surface.  In  the  case  of 
metals,  the  only  value  of  T(k)  of  importance  is  that  determined 
on  the  Fermi  surface. 

c .  The  Diffusion  Equation  for  Electrons  and  Holes 

The  assumed  form  o'"  Boltzmann's  equation  leads  directly  to 
the  fundamental  diffusion  equation  for  charged  particles  in  a 
semiconductor,  which  is  the  basis  for  all  the  analysis  to  date  of 
semiconductor  devices.  If  the  definition  (44)  is  used  in  (42), 
one  finds  the  diffusion  equation  in  the  form 


(48) 


from  which  one  obtains  the  continuity  equation  for  charged  car¬ 
riers  in  the  material.  To  obtain  the  continuity  equation  for 
electrons  in  the  conduction  band,  one  multiplies  by  the  density 
of  states  for  energy  levels  in  the  conduction  band  and  integrates 
over  all  energy  states.  After  proper  manipulation  of  the  terms, 
one  finds 


dn 

at 


(49) 


where  Id  is  a  diffusion  current  density  resulting  from  a  density 
gradient  of  electrons  and  I0  is  a  conduction  current  density 
caused  by  an  electric  field. 

The  current  density  I^  is  obtained  from  the  middle  term  of 
the  right-hand  side  of  (48)  and  can  be  shown  to  have  the  form 
given  by  reference  26 

Id  =  +q  D  vn  (50) 
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where 

q  =  the  charge  on  the  particle 

D  =  an  appropriately  defined  diffusion  constant 

Equation  (50)  is  the  well-known  Fick's  Law  approximation  to  the 
current,  resulting  from  a  density  gradient  of  particles.  Deri¬ 
vations  of  Fick's  Law  exist  in  many  books  and  will  not  be  pre¬ 
sented  here.  Fick's  Law  is  an  approximation  which  depends,  as 
does  the  relaxation  time  assumption,  on  the  particles  not  being 
limited  in  their  time  and  space  movements.  Such  an  assumption 
will  not  be  strictly  applicable  to  physical  sizes  that  are  of 
the  order  of  a  diffusion  length. 

The  current  density  Ig  is  obtained  from  the  last  term  on 
the  right-hand  side  of  (48)  as  is  shown  in  reference  27  and  re¬ 
duces  to  the  form 

!e  =  nq  un  E  (51) 


where 

|in  =  electron  mobility 
E  =  electric  field  intensity 

The  sum  of  equations  (50)  and  (51)  is  the  total  electron 
current  density  In.  Thus  equation  (49)  can  be  put  into  the 
form 


dn 

*3t  “  " 


(52) 


where  the  last  two  terms  in  (49)  have  been  combined  into  In»  A 
similar  expression  can  be  derived  for  the  holes  in  the  valence 
band.  Equation  (52)  for  electrons  and  the  similar  one  for  holes 
form  the  starting  point  for  most  semiconductor  device  analyses. 
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d .  Limitations  of  Diffusion  Theory 

In  the  above  analysis  there  have  been  two  fundamental  as¬ 
sumptions  that  limit  the  applicability  of  this  approach  to  the 
transport  problem.  The  first  is  that  the  definition  of  a  relax¬ 
ation  process  has  been  possible  only  because  of  the  limitation 
to  elastic  or  "almost  elastic"  scattering  processes.  The  sec¬ 
ond  is  that  all  the  above  analyses  have  been  based  on  bulk  prop¬ 
erties  of  the  material.  Surface  and  size  effects  of  the  sample 
were  completely  ignored.  This  latter  fact  may  not  have  been 
stated  explicitly,  but  it  is  implicit  in  the  idea  that  one  may 
define  transport  coefficients  which  are  characteristic  properties 
of  the  material,  and  are  independent  of  the  shape  and  size  of 
the  specimens . 

The  latter  restriction  will  not  be  considered.  Implicit 
in  the  development  of  equations  (30),  (33)/  and  the  definition 
of  a  valid  relaxation  time  is  the  assumption  that  the  particles 
will  have  the  time  and  space  in  which  the  relaxation  process 
may  proceed  uninhibited.  Thus,  the  resulting  diffusion  equa¬ 
tion  (48)  is  a  long-time,  large-scale  approximation  to  the  solu¬ 
tion  of  the  small-scale  problem  of  random  motion  in  the  material. 
In  order  that  the  approximation  should  be  valid,  it  is  neces¬ 
sary  that  one  be  concerned  with  distances  covering  many  mean 
free  paths  and  with  intervals  containing  many  relaxation  times. 

4 .  Remarks 

In  this  discussion  the  fundamental  Boltzmann  equation  (34), 
describing  the  transport  process  in  semiconductors  and  metals 
has  been  discussed.  This  equation  is  valid  for  many  analyses 
but  also,  in  general,  it  is  impossible  to  solve  without  some 
simplifying  assumptions.  The  definition  of  a  relaxation  time 
allows  one  to  solve  the  Boltzmann  equation  and  to  calculate 
steady-state  transport  parameters.  The  calculation  of  a  relax¬ 
ation  time  requires  that  the  scattering  of  particles  be  elastic 
and  that  the  scattering  process  be  independent  of  the  physical 
size  of  the  specimen.  These  two  conditions  are  satisfied  for 
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most  practical  cases.  The  procedure  for  obtaining  the  charge 
continuity  equation  from  the  Boltzmann  equation  was  indicated. 
This  is  one  of  the  fundamental  equations  used  in  semiconductor 
analysis.  Mention  was  also  made  of  the  approximation  made  in 
the  use  of  Fick's  Law  for  diffusion  current.  The  net  purpose 
of  this  discussion  is  to  lay  the  foundation  to  the  diffusion 
theory  approach  to  transport  phenomena  and  to  point  out  the  fact 
that  diffusion  is  a  bulk  phenomenon  and  is  not  applicable  to 
small  regions.  Thus,  if  physical  dimensions  of  specimens  in  the 
direction  of  transport  are  not  greater  than  a  few  diffusion 
lengths,  diffusion  theory  will  not,  without  modification,  de¬ 
scribe  the  transport  process  . 

5 .  Introduction  to  Diffusion  Theory  and  Semiconductor 

Analysis 

Diffusion  theory  has  been  the  basis  of  semiconductor  device 
analysis  since  the  advent  of  such  devices.  Such  analysis,  as 
was  pointed  out  in  the  preceding  discussion,  is  based  on  the 
fact  that  physical  dimensions  of  the  device  are  such  as  to  have 
no  effect  on  the  scattering  process.  In  this  subsection  the  re¬ 
sults  of  applying  diffusion  theory  to  the  study  of  semiconductor 
devices  will  be  discussed.  The  fundamental  equations  for  semi¬ 
conductor  material  analysis  are  written  down.  A  brief  discus¬ 
sion  of  a  digital  computer  code  developed  at  General  Atomics, 
used  to  obtain  the  operating  characteristics  of  a  semiconductor 
diode  in  a  radiation  environment,  is  given.  Finally,  a  discus¬ 
sion  of  general  transistor  models  and,  in  particular,  a  compar¬ 
ison  of  the  three  most  popular  models  after  the  fashion  of 
Hamilton,  et  a  1., concludes  the  subsection. 

6 .  The  Fundamental  Semiconductor  Equations 

Almost  all  device  analyses  have  been  carried  out  in  one 
spatial  dimension;  i.e.,  the  assumed  direction  of  current  flow 
in  the  device.  The  fundamental  equations  used  in  the  analysis 
are : 
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Charge  continuity  equations 
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Total  current  equations 


I  =  qn  u  E  +  qD  (electrons) 

«  ^  V  ^  n  dx  '  ' 
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Poisson's  equation 
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a  ^  _ 


dx 


-7WithEx  =  -f£ 


Charge  neutrality  equation 


p  +  N ,  =  N  +  n 
da 


The  notation  in  equations  (54)  to  (58)  is: 

n  =  electron  density 

nQ  =  equilibrium  electron  density 

p  =  hole  density 

pQ  =  equilibrium  hole  density 

T  =  relaxation  time  for  electrons 
n 

T  =  relaxation  time  for  holes 
P 

D  =  diffusion  constant  for  electrons 
n 

D  =  diffusion  constant  for  holes 
P 

M  =  lectron  mobility 
n 


(53) 

(54) 


(55) 

(56) 


(57) 


(58) 
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Pp  =  hole  mobility 

E  =  electric  field  in  x-direction 

I  =  electron  current  in  x-direction 
n 

I  =  hole  current  in  x-direction 
P 

q  =  electronic  charge 

ijj  t  x  =  electrostatic  potential  at  x 

P  =  charge  density 

€  =  dielectric  constant 

=  density  of  ionized  donor  atoms 
N  =  density  of  ionized  acceptor  atoms 

Equations  (54)  to  (58)  with  the  appropriate  boundary  conditions, 
derived  from  basic  semiconductor  physics,  form  the  starting 
point  for  all  device  analyses. 

The  diffusion  equations  can  be  solved  explicitly  for  semi¬ 
conductor  devices  under  assumed  boundary  conditions,  but  for  any 
device  consisting  of  more  than  one  p-n  junction  the  solution 
becomes  highly  unwieldy.  Therefore,  there  have  been  various 
attempts  to  develop  useful  device  models  that  are  accurate  and 
yet  less  cumbersome.  Before  discussing  some  of  the  better  known 
models,  an  attempt  to  investigate  the  decay  of  excess  charges 
in  a  junction  diode  numerically  will  be  reviewed. 


Digital  Computer  Analysis  of  the  Continuity  Equation 


The  investigation  was  made  under  an  Air  Force  Contract 
AF  29(601)637^  and  reported  in  a  technical  report  AFWL  TR-65-44 
(reference  28).  The  excess  charge  in  this  study  was  generated 
by  an  ionizing  pulse  of  radiation.  The  diffusion  equations  are 
written  in  the  form 
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and  Poisson's  equation 
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(61) 


for  each  side  of  a  p-n  junction* 
the  form 


The  recombination  term  R. 

J 


has 


pn  -  n 
_  __  on  on  / 

Rn - rn -  (n  re<3lon)> 

n  on 


pn  -  p  n 

RP  =  - — -2E-2R  (p  region) 

P  P*op 


(62) 


where 


K  =  the  dielectric  constant 

ANj  =  excess  of  donors  to  acceptors  in  the  j  region 
pQ  =  equilibrium  hole  concentration 
nQ  =  equilibrium  electron  concentration 

The  source  term  was  taken  to  be  proportional  to  the  total 
radiation  received.  The  differential  equations  were  approxi¬ 
mated  by  finite  difference  equations  for  which,  after  some  dif¬ 
ficulty,  solutions  were  obtained.  Some  of  the  difficulties  en¬ 
countered  were  due  to  the  large  range  of  variables,  errors  in¬ 
troduced  when  using  electric  field  and  carrier  density  evaluated 
at  the  same  point,  and  sensitivity  of  solutions  to  minor  changes 
in  boundary  conditions.  Most  of  the  difficulties  were  overcome 
to  a  point  where  the  t ime -dependent  diffusion  equations  were 
solved  for  carrier  densities  and  the  electric  field  in  the 
semiconductor  diodes. 
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Specific  calculations  were  made  for  the  response  of  the 
above  variables  to  a  short  (=  .1  usee)  pulse  of  ionizing  radia¬ 
tions.  These  results  were  then  compared  to  experimental  data  on 
two  particular  diodes.  Fairly  close  agreement  was  found  for 
low  values  of  excess  carrier  densities  with  the  deviations 
greater  at  higher  densities.  It  is  interesting  that,  after  an 
initial  adjustment,  computer  time  steps  greater  than  the  relax¬ 
ation  time  for  the  material  gave  results  agreeing  with  experi¬ 
mental  data. 

Although  the  computer  calculations  reported  in  reference  28 
were  carried  out  only  for  diodes,  the  authors  propose  to  extend 
the  solutions  to  multiregion  devices  in  the  future  and  to  incor¬ 
porate  provisions  for  time-dependent  recombination  rates,  elec¬ 
tric  field  dependent  mobilities,  and  other  phenomena.  This  pro¬ 
cedure  has  produced  acceptable  results  for  transient  analysis  of 
semiconductor  devices.  There  are  numerous  studies  reported  of 
the  transient  response  of  semiconductor  devices  to  ionizing 
radiation  (see  reference  29  and  its  references).  These  are  not 
considered  here . 

7 •  General  Transistor  Models 

It  is  to  be  noted  that  the  set  of  equations  (54)  to  (58) 
subject  to  the  boundary  conditions  at  a  p-n  junction  can  be 
solved  explicitly  for  the  currents  in  a  diode  or  transistor. 

The  diffusion  theory  approach  to  transistor  analysis  is  pre¬ 
sented  in  Chapter  7  of  Middlebrook  (reference  J>1)  or  Chapter  4 
of  Gartner  (reference  50).  These  procedures  develop  the  "law 
of  the  junction"  which  relates  the  minority  carrier  densities 
on  both  sides  of  the  junction  to  the  equilibrium  carrier  densi¬ 
ties  and  the  applied  potential.  The  result  is  that 


P  =  pn  exp  Sf 


(n-material ) 


(65) 


n 


I 


r 


n 


qV 

n  exp  r— • 
p  r  kT 


( p-materia 1 ) 


t64) 


where 

Pn  =  equilibrium  density  of  holes 

=  equilibrium  density  of  electrons 

V  =  applied  voltage  plus  contact  potential  (contact 
potential  is  often  small  enough  to  neglect) 

q  =  electronic  charge  magnitude 

k  =  Boltzmann  constant 

T  =  absolute  temperature 

This  "law  of  the  junction"  provides  the  necessary  boundary  con¬ 
ditions  for  the  solution  of  the  diffusion  equation  in  the  base 
region.  By  assuming  only  diffusion  of  minority  carriers  in  a 
space-charge -neutral  base  region,  the  fundamental  equations  re¬ 
duce  to 


for  the  n-material  base  of  a  p-n-p  transistor.  Solutions  for 
equations  (65)  -  (68)  subject  to  the  law  of  the  junction  at  the 
collector  and  emitter  edges  of  the  base  region  are  given  in 
Middlebrook  for  the  steady-state  and  time-varying  cases  for 
small  signals.  These  solutions  provide  the  emitter  and  col¬ 
lector  currents  under  specified  operating  conditions.  For  the 
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For  this 


exact  form  the  reader  is  referred  to  reference  31** 
discussion,  it  is  sufficient  to  note  .that  the  collector  and 
emitter  currents  can  be  put  in  the  form 


Ie(s)  =  Au(s)  P.(-)  +  Al2(s)  Pc(s) 

(emitter ) 

(69) 

Ic(s)  =  A21(s)  Pe(s)  +  A22(s)  Pc(s) 

(collector ) 

(70) 

where 

Pg  =  hole  density  at  emitter  edge 
Pc  =  hole  density  at  collector  edge 

The  Ajj  are  determined  from  the  solution  of  the  diffusion  equa¬ 
tion  for  holes  subject  to  the  boundary  conditions  given  by 


qV 

(at  collector  junction) 

(71) 

P  =  Pn  exP  kT~ 

qVe 

(at  emitter  junction) 

(72) 

P  =  pn  GXP  kT- 

where 

p  =  equilibrium  density  of  holes  in  base 

Vc  =  collector  to  ‘  >»*e  voltage 

V0  =  emitter  to  base  visit  age 

The  forms  for  Ie  and  I  lead  one  to  treat  the  base  region 
as  a  two-port  network  with  parameters  of  emitter  and  collector 
currents  and  hole  density  at  the  emitter  and  collector 
boundaries.  If  the  area  and  carrier  lifetimes  are  uniform 
throughout  the  base,  then  the  network  will  be  symmetrical,  with 
the  results  that  A^  =  A22  and  A2^  =  A^  •  For  such  a  case,  the 
coefficients  are  given  by 

*See  pages  138  and  140  for  the  steady-state  and  pages  142 
and  i42  for  the  time-varying  case. 
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(75) 


(74) 


where 


W 

A. 


J5 


_T 

P  P 


base  region  width 
base  area 


It  is  apparent  from  (71 )  and  (72)  that  if  Vc  is  large  and  nega¬ 
tive,  then  p  2:  0,  and  if  V  is  large  and  negative  p  ~  o. 

Thus,  by  making  forward  and  inverse  measurements  under  the  above 
conditions,  one  can  define  forward  and  inverse  short-circuit  cur 
rent  ratios  as 


“n  = 


Ic(s) 

iTFT 


and 


pc=0 


Ie(s) 

TJJ7 


Pe=0 


(75) 


In  the  symmetric  case  e*n  =  ,  but  in  practice  the  base  region 

is  seldom  symmetrical.  Consequently,  one  will  obtain,  in  gen- 

eral,  different  values  for  a  and  a.  .  The  measurements  of  a 

n  l  n 

and  are  made  from  direct  current  to  very  high  frequencies, 
and  cun  and  av  are  then  defined  as  the  frequencies  at  which  the 
values  of  «n  and  are  three  decibels  below  the  low-frequency 
values . 

From  reference  50,  one  finds  the  result  of  applying  diffu¬ 
sion  theory  to  the  problem  of  determining  the  time -dependent 
current  in  a  transistor  operated  in  the  common  base  configura¬ 
tion.  It  is 
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where  Iq  is  the  final  value  of  the  collector  current.  Although 
this  result  can  be  used,  it  is  cumbersome.  Therefore,  various 
attempts  have  been  made  to  produce  models  that  are  easier  to  use 
and  give  results  that  are  fairly  accurate.  Some  of  those  models 
will  now  be  discussed,  following  the  review  article  by  Hamilton, 
et  a  1.  (reference  32). 

a .  The  Ebers-Moll  Model 

The  c^'s  and  a^'s  determine  transistor  operating  character¬ 
istics  and  are  easily  measured  quantities;  but  they  have  some¬ 
what  less  than  a  desirable  analytical  form:  e.g.,  for  the 
symmetric  case 


W/L 


a  =  a.  =  sech 
n  l 


(1  +  STp) 
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The  Ebers-Moll  model  assumes  a  more  familiar  form  for  current 
ratios  in  the  forward  and  inverse  regions.  This  model  approxi¬ 
mates  the  a's  by  a  single-pole  function  of  the  complex  frequency 
variable.  In  the  forward  active  regions  this  is  defined  as 
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and  in  the  inverse  active  regions  as 
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This  approximation  allows  each  current  in  the  transistor  to  be 
written  as  the  sum  of  forward  and  reverse  components.  Then  (78) 
and  (79)  yield  equations  for  the  transistor  saturation  currents 
in  the  form 


ie(s)  -  W8) 


(80) 


a^d 


(81) 


The  characteristic  equation  for  the  transistor  in  saturation 
has  the  form 

2 

s  +  (o>n  +  0).)s  +  ayo^l  -  an«.)  =  0. 

b.  The  Charge  Control  Model 

The  essence  of  the  charge  control  model  is  the  relationship 
it  draws  between  terminal  currents  and  total  stored  minority 
charge.  The  procedure  is  to  remove  the  functional  dependence 
on  spatial  corrdinates  by  integrating  the  continuity  equations 
for  the  minority  carrier  charge  density  over  the  volume  of  the 
base  regions.  The  validity  of  such  an  approach  depends  cru¬ 
cially  on  the  steady-state  ratios  of  currents  to  stored  charge 
remaining  the  same  for  all  time.  Performing  the  integration 
leads  to  the  equation 


(82) 
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(85) 


dqBF 

“dt  _iBF  +  qBF  WBF 


where,  by  definition,  the  charge  control  parameter  is 


and  where  I  and  Q  are  steady-state  current  and  charge,  and  i  and 
q  are  time -varying  current  and  charge. 

One  can  define  charge  control  parameters  between  each  cur¬ 
rent,  collector,  emitter  and  base,  and  the  total  charge  for  both 
normal  and  inverse  operation.  The  operations  of  the  transistor 
can  then  be  described  in  terms  of  these  parameters.  Eventually 
the  charge  control  parameter  must  be  related  to  the  <*_,  an,  wn, 
and  Wj  since  these  are  the  measured  quantities.  If  this  rela¬ 
tion  is  made  (see  formulas)  then  it  is  found  that  the  model 
gives  the  same  characteristic  equation  (82)  as  does  the  Ebers- 
Moll  model.  The  charge  control  model  has  perhaps  been  used  more 
than  others  in  the  study  of  radiation  effects.  Under  radiation 
the  total  stored  charge  in  the  base  includes  the  induced  charge. 
Thus  the  effect  of  the  radiation  can  more  easily  be  determined 
since  the  induced  charge  is  more  directly  related  to  the  tran¬ 
sistor  operation. 

c .  The  Linvill  Model 

The  Linvill,  or  lumped  model,  is  a  circuit,  or  network 
analysis,  approach  to  the  problem  of  modelling  semiconductor 
devices.  The  space  dependence  of  the  fundamental  partial  dif¬ 
ferential  equations  is  removed  by  a  group  of  finite  difference 
equations  written  for  a  number  of  points  throughout  the  base 
regions.  Doing  this  gives  the  continuity  equation  at  each  point 
k  in  the  form 

dpk  _  n  [pk+l  ~  2pk  +  pk-l]  (pk  -  pn> 

dt  p  Vhk+1 +  hk)/2  Tp 
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To  establish  circuit -like  elements  to  represent  the  diffusion, 
recombination,  and  storage  process,  Linvill  defines  diff usance 
as 


eAcL 


I 


the  storance  by 


eA 


and  combinance  as 


(85) 
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Then  equation  (84)  is  reduced  to 
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Hckpk  +  Hckpn 


and  gives  the  total  base  current  as 


1b  k+i  CHckPn " 


dpk, 

Hckpk  "  Sh  ^dt  ^ 


(88) 


(89) 


where  N  =  number  of  points  in  the  base  regions. 

The  above  elements  are  given  physical  representation  and 
network  configurations,  then  are  used  to  analyze  the  transistor 
actions.  The  most  popular  configuration  is  a  single  ir  section 
of  the  lumped  elements.  When  a  transistor  is  analyzed  this  way, 
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the  result  is  the  same  characteristic  equation  (82)  as  predicted 

by  the  Ebers-Moll  and  charge  control  models  when  the  lumped 

parameters  are  related  to  the  ai  ,  cd  ,  a.  and  a  . 
r  n  e  i  n 

Thus  each  of  the  above  models  produces  the  same  measure  of 
approximations  to  the  transistor  analysis.  In  all  models  the 
result  is  the  same:  the  space  dependence  has  been  removed  from 
the  fundamental  equations  leaving  a  set  of  ordinary  differential 
equations.  Also  in  each  model,  in  one  way  or  another,  the  steady- 
state  characteristics  are  extrapolated  to  obtain  the  transient 
behavior. 

Based  on  the  diffusion  theory  approximation  of  the  funda¬ 
mental  equations  used  in  semiconductor  analysis,  the  above 
models  are  at  least  as  good  as  any  others  that  have  been  pro¬ 
posed.  Now  the  question  to  be  raised  is:  Is  there  another 
basic  approach  to  the  problem  of  charge  transport  in  semicon¬ 
ductors  that  will  lead  to  equivalent  results  in  devices  where 
diffusion  theory  is  applicable  and  also  valid  in  devices  where 
diffusion  theory  is  not  a  good  approximation?  The  approach 
should  be  correct  in  thin-film  analysis  where  dimensions  in  the 
direction  of  charge  transport  are  of  the  order  of  a  diffusion 
length.  Such  an  approach  has  been  tried  in  the  field  of  semi¬ 
conductors  by  McKelvey  in  a  specific  application  of  the  general 
method  of  invariant  imbedding.  This  approach  will  be  introduced  in 
the  next  subsection. 

8.  Subsection  Summary 

In  this  subsection,  the  diffusion  theory  relations  developed 
in  the  first  subsection  have  been  used  to  see  how  their  applica¬ 
tion  to  semiconductor  devices  has  led  to  the  development  of 
useful  models  for  the  analysis  of  such  devices.  First  a  specific 
computer  attempt  to  solve  the  continuity  equation  under  the  ef¬ 
fect  of  a  short  ionizing  pulse  of  radiation  was  discussed. 

There  was  some  difficulty  in  developing  a  working  program,  but 
finally  calculations  were  made  for  specific  diodes  which  agreed 
with  experimental  data.  Next,  models  developed  for  the  general 
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analysis  of  transistor  operation  (including  the  Ebers-Moll, 
charge  control,  and  lumped  models)  were  reviewed. 

The  most  important  point  here  is  not  that  each  of  the 
models  is  just  an  approximation  to  the  fundamental  diffusion 
equations  but  that  it  is  impossible  to  extend  these  models  into 
thin-film  semiconductor  devices  where  diffusion  theory  on  which 
they  rest  does  not  accurately  describe  the  physical  process. 

An  attempt  to  do  this  through  the  general  approach  of  invariant 
imbedding  is  the  purpose  of  the  rest  of  this  section. 

9*  Introduction  to  the  Application  of  Invariant  Imbedding 

Invariant  imbedding  essentially  reverses  the  classical  ap¬ 
proach  to  solving  transport  problems.  Historically,  the  solu¬ 
tions  of  a  particular  transport  problem  are  obtained  by  consid¬ 
ering  special  limiting  cases  of  the  general  formulation  when 
necessary  restricting  and  simplifying  approximations  are  made. 
This  leads  to  a  system  of  linear  equations  which  one  then  solves. 
A  good  example  of  such  a  process  is  the  diffusion  theory  applica¬ 
tion  to  semiconductor  device  operation.  The  invariant  imbedding 
approach  is  to  state  the  particular  problem,  then  to  phrase  it 
in  a  more  general  manner  and  apply  mathematical  invariance  prin¬ 
ciples  in  a  systematic  fashion  to  reduce  the  problem  t^  the 
iteration  of  nonlinear  transformations.  This  procedure  then 
leads  to  new  computational  forms  that  seem  well  fitted  to  the 
capabilities  of  digital  computers.  It  is  not  the  desire  of  this 
study  to  formulate  the  general  procedure  but  merely  to  apply  it 
to  a  specific  problem.  In  this  vein  of  thinking  the  work  of 
McKelvey  will  first  be  reviewed.  He  has  applied  invariant  im¬ 
bedding  to  semiconductor  analysis.  Next,  the  general  time- 
dependent  problem  will  be  formulated  and  the  general  procedure 
for  obtaining  a  solution  will  be  indicated. 

10.  Time  Independent  Reflection  and  Transmission  Coefficients 

Invariant  imbedding  started  as  essentially  a  particle  count¬ 
ing  technique.  This  technique  is  used  by  McKelvey  in  his  earlier 
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works  (references  33 »  3^/  and  35)  to  study  recombination  rates 
of  injected  carriers  in  semiconductors  due  10  volume,  and  at 
dislocation  edges  in  the  material.  In  these  papers  the  partic?^e 
counting  method  was  used  to  derive  reflection  and  transmission 
coefficients  for  the  surfaces,  volumes,  or  dislocations  consid¬ 
ered.  These  coefficients  were  then  used  to  determine  appropri¬ 
ate  boundary  conditions  for  the  solution  of  the  diffusion  equa¬ 
tions  in  the  material.  In  these  earlier  papers  the  reflection 
and  transmission  coefficients  were  related  to  previously  defined 
physical  properties  and  are  not  developed  in  the  manner  of  in¬ 
variant  imbedding.  Thus,  these  articles,  although  providing 
background  for  McKe Ivey's  work,  are  not  essential  to  an  under¬ 
standing  of  his  application  of  invariant  imbedding. 

For  purposes  of  introduction  the  formulation  for  the  funda¬ 
mental  differential  equations  for  the  reflection  and  transmis¬ 
sion  coefficients  as  developed  by  McKeivey  in  reference  36  are 
given  below.  He  considers  the  one-dimensional  flow  of  particles 
through  a  given  medium  as  presented  in  figure  . 


Figure  j4.  Th^  Flukes  4  n  «*  Stratified  Medium 


To  develop  the  equtet j-c.-nm  a  **U1  of  length  x  +  dx  is 
considered;  reflection  afrd  tranun,^  ,:on  coefficients  R(x)  and 
T(x)  are  defined  as  functions  of  *.ength  x.  Then  the  reflec¬ 

tion  and  transmission  f*  *  *  ckness  dx  is  proportional  to 
dx  and  one  writes 
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R ( dx )  =  kdx 

(90) 

T(dx)  =  1  -  ft  dx 

(91) 

The  parameter  ft  includes  both  scattering  and  absorption;  there¬ 
fore  ft  >  k,  the  equality  holding  when  there  is  no  true  absorp¬ 
tion  but  only  scattering.  A  is  the  source  term  of  incident 
particle  flux,  and  the  are  fluxes  across  the  boundaries  in 
the  direction  indicated.  One  can  obtain  the  equation  by  one  of 
two  equivalent  methods:  the  first  by  considering  multiple  re¬ 
flections  and  transmissions  as  was  originally  done  in  reference 
37  or,  second,  by  bulk  considerations.  The  latter  method  will 
be  used.  From  a  consideration  of  figure  34  the  following  equa¬ 


tion  can  be  written 

F^  =  A(l  -  ft  dx)  +  kF^  dx  (92) 

=  R(x)P1  (9?) 

F^  =  Ak  dx  +  (l  -  ft  dx)F^  (94) 

F,  =  t(x)fl  (93) 


From  equation  (93)#  by  straightforward  procedures  one  can 
obtain 
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aMftt 


J-  -flWfcsJjSS 


AT(x ) 


l  -  g  dx 
1  -  R(x)  k  dx 


(99) 


Since  is  the  reflected  flux  at  (x  +  dx),  the  reflection  co¬ 
efficient,  R(x  +  dx),  is,  from  (97) 

R(x  +  dx)  =  .  k  a*  +  (100) 

Then,  neglecting  second-order  terms,  one  obtains  for  ^100) 


dR  ,  R(x  +.  dx),  -  Rjxl  =  k  .  2or(x)  +  kR(x)2 


as 


(101) 


which  provides  a  form  for  the  reflection  coefficient  as 

V 

—  sinh  qx 

cosh  qx  +  |  Sinh  qx 
2  2 

where  q  =  (<y  -  k  ) 

In  a  similar  manner  one  can  obtain  the  equation  for  T  from 

t(x  +  dx)  -  -  T(x)  R(x)d^)dx  (10?) 

leading  to  the  differential  equation 

£  =  T(X  +  ^  ^  =  T(x)[R(x)k  -  «]  (104) 

Therefore,  although  the  equation  ( 101 )  for  R  involved  only  R 
itself,  T  cannot  be  found  from  (105)  without  knowing  R.  If  R 
from  (102)  is  used  one  can  eventually  show  that  T(x)  is  given 
by 


R(x) 


98 


I 


T(x)  = 


cosh  (qx)  +  —  sinh  qx 


(105) 


a .  Analysis  of  a  p-n  Photodiode 

McKelvey,  after  arriving  at  (102)  and  (105),  relates  k  and 
Ot  to  the  usual  semiconductor  parameters  of  mean  free  path,  mean 
carrier  lifetime,  and  mean  thermal  velocity.  This  will  not  be 
repeated  here.  McKelvey  then  goes  on  to  use  R(x)  and  T(x)  as 
found  to  analyze  the  operation  of  a  p-n  photodiode  on  the  flux 
basis.  The  expression  for  the  differential  flux  dF  across  the 
junction  due  to  a  differential  source  term  dA  given  by 


dF (y  )  =  $(x,  y )dA 


(106) 


where  #(x,  y)  takes  on  the  identity  of  a  Green's  function  for 
the  problem  and  y  is  a  parameter  of  the  material  defined  by 


2  \  _  m  2  /  scattering  mean  free  path  \ 
T  L  Y  T  '■diffusion  length  of  the  particles' 


(107) 


For  the  diode  a  bulk  generation  rate  of  carriers  g(x)  is  assumed 
Thus  dA  =  g(x)  dx,  and  in  order  to  find  the  total  flux  across 
the  junction  one  integrates  (106),  giving 


VA 

F(y)  =  ^  g(x)*(x,  Y  )dx 


(108) 


where  the  surface  is  at  x  =  0,  and  the  junction  is  at  x  =  d. 
The  function  $  has  the  form 


$(x,  Y)  = 


[ 


T  ( x )  Ro 


R(x)ir  1  ,  R(x ) 

T(x )  lt(x  )  T(x ) 


T(a )  [tT3T7  ro  T 


-  R. 


®  " 


1  +  R, 


R(a  -  x 


o  T(a  -  x 


(109) 
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with  R  =  surface  reflection  coefficient,  and  "a"  the  slab  width, 
o 

This  equation  can  be  put  into  other  forms  using  (105)  and  ( 102 ) 
where  y  is  introduced  through  k  and  a. 

The  integral  equation  (108)  is  the  fundamental  equation 
involved  in  the  analysis  of  the  problem  using  the  flux  or  in¬ 
variant  imbedding  method.  McKelvey  analyzed  the  photodiode  for 
the  source  term 

g(x)  =  90*(x  -  xQ)  0  *  x  *  a  (110) 

where  fi(x)  is  the  Dirac  6  function,  and  showed  that  the  results 
of  such  an  analysis  agree  with  diffusion  theory  results  of  the 
problem  in  the  limit  where  the  invariant  imbedding  method  and 
diffusion  theory  are  the  same.  This  limit  occurs  when  y  «  1, 
as  can  be  seen  from  the  definition 


Y 


_  2  X 

‘  J  L 


(111) 


2  A 

If  y  —  «  1,  then  the  mean  free  path  is  short  compared  to  the 
diffusion  length  L.  Or,  if  L  =  cr  where  c  =  average  thermal 
velocity  and  r  =  mean  carrier  lifetime,  the  condition  implies 
a  combination  of  short  mean  free  path  or  long  lifetime  or  high 
velocity  that  will  make  y  «  1.  Then  McKelvey  established  in 
this  case  that  invariant  imbedding  and  diffusion  theory  agree 
if  diffusion  theory  is  applicable.  It  was  shown  in  the  first 
subsection  that  diffusion  theory  was  applicable  if  the  physical 
size  were  large  compared  to  a  mean  free  path.  Therefore,  if, 
for  example,  the  surface  layer  of  the  material  were  made  very 
thin,  the  diffusion  theory  analysis  would  be  invalid  but  invar¬ 
iant  imbedding  could  still  be  used  to  analyze  the  operations  of 
the  junctions.  In  a  following  paper  (reference  38)  McKelvey 
continued  the  invariant  imbedding  analysis  of  a  p-n  junction. 
Alternative  forms  for  T  and  R  are  defined  which  simplify  the 
arithmetic  somewhat.  An  expression  for  the  total  current  across 
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a  p-n  junction  is  derived  using  invariant  imbedding.  The  cur¬ 
rent  is  given  by 


I  =  q 


DR  DR  l^2 

np<-2-=&)  ♦  P„(-|-=E) 

*  n  p 


(“1>S 


- 1) 


(112) 


where 

R  =  the  reflection  coefficient  for  n-material  of  infinite 
°°n  length 

R  =  the  reflection  coefficient  for  p-materials  of  infinite 

cor}  • 

*  length 

and  the  other  quantities  are  as  defined  in  division  5»  With 

1/2  1/2 

the  exception  of  (R^)  and  (R  )  this  equation  is  identi¬ 
cal  to  the  current  found  from  diffusion  theory  given  by 


I  =  q 


D  D 

n  _£  +  p  _E 

P  L  pn  L 
n  p 


(exp  -  !) 


(115) 


where 


L  =  \/  T  D 
P  P  P 


and 


Ln  = 


The  R_'s  are  given  by  the  equation 


R  =  1  +  2y2  -  2y(l  +  Y2) 


1/2 


(114) 


with  y  <<  1 #  the  R^  —  1  and  invariant  imbedding  agrees  with 
diffusion  theory  as  previously  discussed.  Where  this  condition 
fails,  R  <  1  and  invariant  imbedding  predict  a  smaller  satura¬ 
tion  current  than  is  given  by  diffusion  theory.  The  article 
then  continues  with  the  invariant  imbedding  approach  to  the 
analysis  of  a  p-n  photovoltaic  cell,  but  that  will  not  be  dis¬ 
cussed  here. 
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McKelvey's  final  article  to  date  (reference  39)  introduces 
the  effect  of  an  electric  field  into  the  above  analysis.  This 
considerably  complicates  the  problem  because  it  becomes  neces¬ 
sary  to  consider  reflection  and  transmission  coefficients  in 
the  direction  of  the  electric  field  and  in  opposition  to  it. 

More  bookkeeping  is  required  in  this  latter  case  but  in  no  way 
alters  the  basic  idea  behind  the  applications  of  invariant 
imbedding . 

McKelvey's  work  has  been  treated  rather  fully  because  it 
is  a  good  example  of  the  application  of  the  invariant  imbedding 
method  to  the  study  of  semiconducting  devices. 

11.  Remarks  on  Time  Dependent  Invariant  Imbedding 

McKelvey's  series  of  papers  do  not  permit  the  inclusion  of 
time-varying  fluxes.  There  exists,  however,  time  dependent 
formulation  of  the  invariant  imbedding  formalism,  which  has  been 
applied  to  a  study  of  neutron  transport  in  a  rod  (references  40 
and  4l).  The  results  of  the  treatment  are  expressions  for  two 
functions,  r(x,  t)  and  T(x,  t),  the  reflection  and  transmission 
coefficients.  The  details  leading  to  these  expressions  will  not 
be  repeated  here  for  two  reasons:  the  results  are  not  strictly 
applicable  to  a  study  of  semiconducting  devices,  and  there  is 
rather  considerable  algebraic  detail.  It  is  found,  however,  in 
the  above-mentioned  references,  that  the  reflection  and  trans¬ 
mission  coefficients  satisfy  equations  of  the  form 

H  +  §  If  =  a(x)B(X|  +  2£f(x'  O  "  l]a(x)R(x,  t)  + 

t 

a(x)B(x,  t)  $  r)  R(Xt  t  _  T)dT  (115) 

o 


and 
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|  =  [F(x,  t)  -  1  ]  a  ( x )  T  ( x  ,  t)  + 

t 

»(*)B(x.  t)$^%^T(x,  t-T)lT  (116) 
o 

In  the  above  relations,  a(x) Ax  is  the  probability  that  an  inter¬ 
action  will  take  place  in  the  medium  under  discussion  in  the 
region  (x,  x  +  Ax);  F(x,  t)  is  the  internal  generated  flux;  the 
factor  unity  is  the  assumed  external  flux  while  B(x,  t)  is  pro¬ 
portional  to  the  flux  generated  due  to  particle  interactions. 

If  both  R  and  T  are  independent  of  time  and  if  F  =  0,  one 
readily  finds  that  the  above  relations  reduce  to 

H  =  ct(x)b(x)  -  2a(x)R(x)  +  a(x)B(x)R2(x)  (117) 

and 

=  T(x)[-c(x)  +  cj(x)r(x)b(x)  ]  (118) 

These  last  two  relations  can  obviously  be  compared  term  by  term 
with  McKelvey's  fundamental  equations.  Indeed,  to  show  that 
such  a  comparison  is  possible  is  the  sole  point  of  the  remarks 
made  in  this  paragraph  to  this  point . 

These  studies  have  progressed  to  the  point  where  a  time  de¬ 
pendent  generalization  of  McKelvey's  work  will  soon  be  available. 
The  internal  fluxes  arising  from  the  interaction  of  an  external 
gamma  pulse  with  the  material  have  been  specified.  The  terms 
in  the  equations  which  describe  the  absorption  of  the  charge 
carriers  due  to  the  existence  of  trapping  states  in  the  material 
and  electron-hole  recombination  need  further  work.  When  these 
equations  are  written  down,  the  following  results  are  expected. 
There  will  exist  a  diffusion  type  of  equation  in  each  of  the 
(thin)  layers  of  a  (microelectronic)  circuit  which  will  contain 
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exralicitly  the  effects  of  layer  width  on  the  total  time -dependent 
current.  When  such  equations  are  available,  generalized  models 
of  thin-film  devices  will  proceed  in  a  straightforward  way,  after 
the  fashion  discussed  in  Hamilton  et  al. (reference  32). 
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SECTION  IV 

GENERAL  CONCLUSIONS  AND  RECOMMENDATIONS 

r 


From  the  theoretical  and  experimental  studies  on  this 
project,  the  following  conclusions  can  be  made: 

(1)  The  most  important  contributing  factors  to  transient 
radiation  effects  in  microcircuits  are  the  active  junction  re¬ 
gions  of  the  transistor  elements.  Since  the  current  pulse  re¬ 
sulting  from  the  radiation  is  directly  proportional  to  the  active 
volume,  it  would  seem  reasonable  to  keep  these  regions  as  small 
as  possible.  This  implies  the  minimization  of  the  emitter  base 
junction  and  the  collector  junction  areas,  depletion  regions, 
and  minority  carrier  lifetimes.  This  further  implies  that  the 
smallest  possible  transistor  devices  consisting  of  the  required 
power  and  current  levels  should  be  used  in  microcircuits. 

Minority  carrier  lifetimes  can  be  reduced  by  gold  doping,  which 
also  implies  the  use  of  high-frequency  devices.  The  depletion 
width  can  be  minimized  by  using  a  low  collector  voltage. 

(2)  The  microcircuit  should  be  designed  to  use  the  lowest 
impedance  consistent  with  other  design  criteria.  This  require¬ 
ment  is,  of  course,  somewhat  contradictory  to  the  previous  rec¬ 
ommendation  of  reducing  the  active  volume  region  for  the  tran¬ 
sistor  elements,  so  some  compromise  between  these  two  factors 
must  be  made  to  minimize  the  radiation  effect. 

(3)  Experimental  data,  with  the  theory  developed  here  for 
use  with  the  flash  X  ray,  offer  a  method  of  obtaining  minority 
carrier  lifetimes,  depletion  layer  width,  junction  area,  diffu¬ 
sion  length,  and  active  volume  when  the  absorbed  dose  and  the 
ionization  efficiency  are  known  in  the  transistor  elements  of 

a  microcircuit.  Conversely,  given  the  above  parameters,  it  is 
possible  to  determine  the  primary  photocurrent  from  the  micro- 
circuit  junction  under  a  variety  of  radiation  conditions. 
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(4)  It  has  been  shown  that  the  fundamental  Boltzmann  equa¬ 
tion  describing  the  charge  transport  in  semiconductors  can  be 
said  to  modify  the  usual  diffusion  or  continuity  equation  to  ob¬ 
tain  more  accurate  circuit  models  of  transistor  elements  in 
microcircuits.  This  will  result  in  changes  in  the  Ebers-Moll 
model,  charge  control  model,  and  the  Linv.ill  model.  Further 
work  will  be  required  to  complete  the  changes  in  these  models, 
but  adaptations  to  thin- layered  microcircuit  devices  has  been 
shown  to  be  feasible . 

(5)  In  general,  the  Motorola  XC-201  is  less  sensitive  than 
the  Motorola  MC-201  to  transient  radiation.  The  XC-201  and 
MC-201  are  four-input  integrated  circuits  "nand"  or  "nor"  gates. 
The  XC-201  has  its  components  completely  isolated  from  each 
other  without  the  use  of  a  back  biased  isolation  junction  diode. 

(6)  Testing  of  the  Fairchild  mL-903  RTL  integrated  circuit 
indicated  that  the  substrate-collector  junction  introduced  a 
diffusion  tail  in  the  transient  radiation  response  which  can  be 
reduced  by  elimination  of  the  diode  isolation  junction.  The 
successful  mocking  up  of  the  mL-905  showed  this  technique  to  be 
a  useful  research  technique. 

(7)  Mock-up  circuits  for  the  TI  SN-513  "nor"  or  "nand" 
integrated  logic  circuits  required  the  inclusion  of  the  parasitic 
elements  for  accurate  simulation. 

(8)  The  TI  SN-511  and  SN-5HA  integrated  circuit  flip- 
flops  were  irradiated.  The  irradiation  tended  to  try  to  switch 
the  flip-flop,  but  it  restored  itself  to  its  original  state. 

(9)  The  Crystalonics  2N3084  field  effect  transistor  was 
tested  on  the  600-KV  flash  X  ray.  The  test  results  indicate 
that  the  transient  peak  of  I_  increases  as  V__  decreases.  The 

(j  bo 

test  results  also  indicate  that  the  transient  peak  of  1^  de¬ 
creases  with  increasing  VDg  around  the  VGg  =  0  volt  level,  and 
that  ID  increases  with  increasing  VDg  at  the  VQg  =  4  volts  level. 

(10)  Diffused  resistors  have  a  large  transient  current 
similar  to  that  observed  in  many  microcircuits.  The  current 
pulse  is  probably  caused  by  the  isolation  junction. 
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(11)  The  transient  radiation  response  of  thin-film  nichrome 
resistors  is  primarily  due  to  air  ionization  effects. 

( 12 )  Thin-film  Hall  devices  were  found  to  be  highly  radia¬ 
tion  resistant  as  was  predicted  by  the  theory.  There  were  no 
radiation  effects  observed  for  dose  rates  up  to  10^  roentgens/ 
sec  when  the  devices  were  operated  in  the  constant  voltage  mode. 
This  operating  mode  is  therefore  preferred  for  any  Hall  device 
where  gamma  radiation  environments  might  be  encountered  .  Tran¬ 
sient  radiation  effects  were  observed  for  other  modes  (  c/ Ra¬ 
tion  as  predicted  by  the  theory. 

(13)  The  Signet ics  SU-315K  integrated  "m  circuit  is  in¬ 
sensitive  to  gamma  radiation  while  in  the  logic  "1"  state--at 
least  up  to  2  x  10°  roentgens/sec.  It  is  sensitive  to  radiation 
at  this  level  while  in  the  logic  M0"  state.  At  high  radiation 
levels  the  SU-jJ15K  becomes  sensitive  to  gamma  radiation  irrespec¬ 
tive  of  its  logic  state. 
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APPENDIX  I 

DERIVATION  OF  p-n  JUNCTION  TRANSIENT  CURRENT 

In  order  to  evaluate  an  expression  for  the  transient  re¬ 
sponse  of  a  p-n  junction  several  assumptions  have  to  be  made. 
These  assumptions  are: 

(1)  The  voltage  applied  across  the  p-n  junction  is  as¬ 
sumed  to  occur  only  across  the  area  in  the  immediate  vicinity 
of  the  junction.  The  field  available  for  drift  within  the  main 
body  of  the  p-region  or  n-region  is  therefore  very  weak. 

(2)  The  diode  has  a  one -dimensional  carrier  flow. 

(3)  The  diode  is  uniformly  doped. 

(4)  The  ionizing  radiation  produces  uniformly  excess 
electron-hole  pairs.  These  excess  electron-hole  pairs  produced 
are  assumed  to  be  only  in  sufficient  amount  to  significantly 
alter  only  the  minority  carriers  in  each  region  as  shown  in 
figure  2  .  When  considering  the  transient  response  of  a  p-n 
junction,  one  runs  into  the  problem  of  hole  injection  into  an 
n-region,  or  electron  injection  into  a  p-region  (reference  6). 
If  holes  injected  into  an  n-region  are  considered,  one  can  ex¬ 
press  in  words  the  conservation  of  holes  at  a  fixed  point  in 
space  as 


time  rate  of 
increase  in 
hole  density 


thermal  generation 
rate  of  holes 


recombination 
rate  of  holes 


outflow  1  + 
of  holes_ 


radiation 
generation 
rate  of  holes 


(119) 


The  first  term  on  the  right-hand  side  of  the  above  equation  is 
o 

equal  to  rn£  and  the  second  term  is  equal  to  rnp.  If  it  is  as¬ 
sumed  that  only  the  minority  carrier  density  is  altered 
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appreciably  in  excess  of  the  equilibrium  densities 
aue  to  hole  injection,  one  can  write 


p„  =  +  6p„ 

n  no  n 


n  =  n 
n  no 


no 


and  n 

no 


(120) 

(121) 


Hence,  if  one  neglects  momentarily  the  radiation  term,  the  first 
two  terms  on  the  right-hand  side  of  equation  (119)  can  be 
written 


d(*pj 


dt 


—  =  rnf  -  rnp  =  rnf  -  r  n  (p  +  ftp  ) 

i  r  l  no v rno  *n ' 


(122) 


Since  rnf  is  equal  to  rn  p  .  equation  (122)  becomes 
i  no  no  v  ' 


3t ( f pn )  =  -rnno*Pn 


(125) 


Solving  for  ftpn  gives 

-rn  t  ~t/  T_ 

*Pn  =  («P)oe  =  («p)Qe  p  (124) 

where  one  defines  rn  _  =  1  /r  ,  and  (ftp)  is  the  value  of  ftp 

no  p  x  c'o  n 

at  t  =  0.  Solving  for  ftpn  out  of  equation  (120)  and  putting  it 
into  equation  (123)  obtains 


d 

It 


(*pn) 


”r  nno^n 


Pno  ^ 


(125) 


The  outflow  of  holes  can  be  found  by  examining  a  cube  of 
widths  Ax,  Ay,  Az .  Net  inflow  into  cube  = 


1 

'_q 


bJv  1  6Jz 

^-1-  Ay  +  ±  ■*-=■  Az 
ay  1  q  8z 


+ 


1  ^x 
q  8x 


(126) 
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Or»  consider ing  a  unit  cube,  one  can  rewrite  equation  (126)  as 


oy  oz 


(127) 


where 


is  the  current  density. 


The  continuity  equation  for  holes  injected  into  an  n-region 
becomes  from  the  above  derivation 


opn(x,  t) 

at 


no 


-  Pn(x. 


t) 


•  Jp  +  g(t) 


(128) 


where  g(t)  is  the  carrier  generation  rate  due  to  radiation. 
Similarly  the  continuity  equation  for  electrons  can  be  written 


Bn  (x,  t) 
~~^t - 


V 


3n  +  g(t)  . 


(129) 


The  change  in  signs  arises  due  to  the  negative  charge  of  elec¬ 
trons.  Since  the  diffusion  component  is  proportional  to  a  car¬ 
rier  density  gradient,  the  diffusion  equations  for  electrons 
and  holes  are 


jp(x.  t) 


dpn(x-  t) 
-’"p-Tx - 


(130) 


j 


n 


t) 


dnp(x,  t) 
dx 


(131) 


The  transient  response  of  a  diode  is  made  up  of  two  com¬ 
ponents:  One  component  is  a  prompt  photocurrent,  i  (t )  ,  which 
is  produced  by  carriers  generated  in  the  junction  transition 
region.  Since  this  component  is  collected  in  a  few  nanoseconds, 
it  can  be  treated  as  having  no  time  delay  with  respect  to  the 
radiation  pulse: 
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ip(t)  =  qA  Wt  g(t) 


(152) 


The  second  component  is  a  delayed  component  of  photocurrent, 
i^(t),  developing  from  the  diffusion  of  carriers  generated  in 
the  vicinity  of  the  junction: 


op 


id(t)  =  qADp^n 


x=0+ 


dn 

qAD  t— 2. 

+  w  MV 


n  ox 


x=0- 


(155) 


If  one  uses  the  boundary  conditions 

Pn(0+-  t)  =  PnoeqvAT,  np(0- ,  t)  =  npoe^AT  (154) 

lim|pn(s,  t)|  <  «,  lim  lnp(x ,  t)|  <  *  (135) 

equations  (128)  and  (129)  can  be  solved  by  use  of  Laplace 
transforms.  By  using  equation  (l^l)  one  can  write  equation 
(129)  as 


dnp(x,  t) 
dt 


a2n 


-  nr,(X-  t) 


+  g(t)  (156) 


n 


Taking  the  Laplace  transform  of  this  equation  one  obtains 


*2  n„(x,  0) 

sNp(x,  s)  -  npo(x,  0)  =  Dn  Np(x,  s)  +  P° - 


N  (x,  s) 

-£-= -  +  G(s) 


(157) 


n 


or 
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d^N  (x,  s)  .  , 

°  =  — 2-2 - 5-(s  +7->Vx'  s>  + 

dx  n  n  r 


1_ 

D 


n 


n  (x ,  o ) 


7 -  +  G(s)  +  npQ(x,  0) 


This  equation  is  of  the  form 


(138) 

i 

/ 


0  =  y  +  Ay  +  B 


(139) 


The  solution  to  this  equation  is 


y  =  Np(x,  s) 


-  G( 


s  + 


!  .  exp  x  /-L(  s  +  i-) 


V 


n  n 


n^fx ,  o) 

+ -E°_ -  (i40) 


Similarly  it  can  be  shown 


Pn(x,  s) 


=  G(s) 


S  +  — 

Tpu 


1  - 


exp  -x  h-(s  +7-) 
v  P  P  J 


P  (x#  o) 
rno v  ' 


(141) 


Upon  substituting  equations  (l40)  and  (l4l)  into  the  Laplace 
transforms  of  equations  (132)  and  (133)  the  equation  for  the 
transient  current  is  found  to  be 


lpp(s)  =  qA  G(S) 


wt  ♦ 


/D  r 
n  n 


'1  +  ST, 


n 


^/D  t 

P  P 

^1  +  ST 


PJ 


(142) 


where  I  (s)  is  radiation-induced  primary  photocurrent  Laplace 
PP 

transforms  composed  of  depletion  and  diffusion  components. 

If  G(s)  is  a  pulse  of  magnitude  g  and  duration  tQ  the  in¬ 
verse  transform  of  Ipp(s)  becomes 
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(144) 


t  >  t 


The  above  results  are  similar  to  those  obtained  by  previous 
investigators  (references  2,  7,  8,  9/  10,  11 ,  13 ,  and  17). 

In  planar  and  mesa  transistors  most  of  the  primary  photo¬ 
current  originates  in  the  collector.  By  making  this  assumption 
for  an  n-p-n  transistor  the  rest  of  the  terms  in  the  above  equa¬ 
tion  may  be  neglected,  which  results  in 


ippC*  )  =  <1*5 


Vd 


T 

P  P 


u(t )  erf  »Jt/T  - 

Xr 

i 


(1^5) 


where 


the  electron  charge  =  1.6  x  10”^  coulombs 

the  base-to-collector  junction  area 

electron-hole  pair  generation  rate 

width  of  the  depletion  region 

the  diffusion  constant  for  holes  in  an  n-type 
material 

the  lifetime  of  holes  in  an  n-type  material 
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D  =  the  diffusion  constant  for  electrons  in  a  p-type 
n  material 


r  =  the  lifetime  of  electrons  in  a  p-type  material 
n  _ 

L  =  ,/d  t  =  the  diffusion  for  electrons  in  a  p-type 
n  n  n 

material 


,/D  t  =  the  diffusion  length  for  holes  in  an  n-type 
P  P 

material 


u(t)  and  u(t  -  t  )  =  unit  step  functions  turned  on  at 

t  =  0  and  t  =  t  ,  respectively. 


Figure  35  shows  a  typical  response  one  would  expect  for 


various  approximations  of  , 


V  and  V 
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Figure  35*  Typical  Transient  Radiation  Response  of  a 

p-n  Diode  for  Various  Values  of  L  ,  L  ,  and  W^. 

p  n  t 
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APPENDIX  II 


DEPENDENCE  OF  DEPLETION  REGION  WIDTH  UPON  VOLTAGE 


For  a  grown  junction  n-p-n  transistor,  the  base  conductiv¬ 
ity  is  usually  at  least  an  order  of  magnitude  larger  than  the 
collector.  This  means  that  the  positively  charged  donor  ion 
concentration  in  the  collector  N^  is  smaller  than  the  negatively 
charged  acceptor  ion  concentration  in  the  base  N  .  The  impurity 
concentration  near  the  junction  may  then  be  approximated  by  a 
charge  density  which  changes  linearly  from  a  predominantly  p- 
doped  region  to  a  predominantly  n-doped  region  as  shown  in  fig¬ 
ure  36a  (references  6  and  12):  that  is, 

Nd  -  Na  =  ax  ( 146) 


In  figure  36,  the  origin  has  been  arbitrarily  chosen  at 
the  point  of  zero  net  impurity  concentration.  When  the  col¬ 
lector  is  biased  negatively  the  mobile  carriers  near  the  junc¬ 
tion  are  swept  away,  and  the  remaining  immobile  impurity  ions 
result  in  a  net  charge  density  in  the  region  of  depletion  of 
the  mobile  carriers.  This  net  charge  density  p  will  give  a 
change  in  electrostatic  potential  V  as  given  by  Poisson's 
equation 


2 

V^v 


(cgs  units) 


(147) 


where  €  is  the  absolute  dielectric  constant  of  the  material. 

In  figure  36c,  the  net  charge  density  has  been  linearly  approx¬ 
imated,  and  it  is  given  by  the  equation 


p  =  +qax 


(148) 
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Figure  36.  Characteristic  Results  of  Charges  Across  the 
Transition  Region  at  the  Base  Collector  Junction 
of  a  Grown  Junction  Transistor 


By  using  equation  (l48)  and  assuming  one-dimensional  analysis, 
equation  ( 147  )  becomes 


(149) 


Integrating  with  respect  to  x  gives 

c 

where  E  is  the  electric  field. 

Using  the  boundary  conditions  that  the  electric  field 
(dv/dx)  goes  to  zero  at  the  edges  of  the  depletion  region 
(±^  t),  the  constant  of  integration  (C^)  must  be 

ci=ff -«t  (151) 

and  therefore 


„  _  dv  _  -47rqa  xZ 

_E  -  -S7 - (  T 


+  c. 


(152) 


Integrating  again  with  respect  to  x  gives 
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_  -2vqa 
€ 
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wi 

5“ 


+  c. 


(155) 


Using  the  boundary  condition  that  the  electrostatic  potential 
is  zero  at  the  origin,  the  constant  of  integration  (C2)  must  be 
zero.  This  gives 
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-27rqa 


x2  W 

T 


2' 

t 

TT 


(154) 
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But  the  total  potential  change  Vfc  across  the  junction  is  given 
by 


(155) 


which  from  equation  (154)  is  found  to  be 


(156) 


From  the  above  equation  it  is  seen  that 


(157) 


The  ionized  charge  in  one  side  of  the  space  charge  volume 
is  obtained  by  integrating  equation  (l48)  from  x  =  0  to 
x  =  Wt/2;  this  gives 


Wt/2 

Q  =  \  qa  x  dx 
o 


qa  W2 

— B- 


(158) 


Using  equation  (157)*  equation  (158)  becomes 


(159) 


The  incremental  capacitance  per  unit  area  ( -Q/oVt )  is  the 
value  required  for  the  small  signal  equivalent  circuit,  which 
from  equation  (159)  is 


C 


nQ  _  ( 

bv^  ~  W 


1/3 


(160) 
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With  the  use  of  equation  (157)/  the  above  equation  may  also  be 
written  as 


(161) 


The  above  capacitance  is  for  a  unit  junction,  and  therefore 
the  capacitance  measured  on  a  device  must  be  divided  by  its 
junction  area.  In  equation  (i60)  it  can  be  seen  that  capacitance 

I  /‘Z 

varies  proportionally  ( l/vt )  .  Depending  on  the  processes  of 

transistor  manufacturing,  doping,  etc.,  the  capacitance  can  vary 

from  (l/v.  )^//^  to  ( l/V.  )  .  From  equation  (157)  it  can  be  seen 

t  t  l/3 

that  the  depletion  width  W,  is  proportional  to  (v.  )  ,  and  it 

1  /f>  l  /?  t 

also  can  vary  from  (vfc)  '  to  (vt)  '  . 

If  L  ,  L  ,  and  W.  are  about  the  same  value,  one  should  be 
p  n  t 

able  to  see  the  transient  radiation  effect  of  the  depletion  com¬ 
ponent  as  shown  in  figure  35*  In  equation  (157)  it  is  shown 

1/3 

that  the  depletion  width  varies  as  V  .  Thus,  as  the  reverse 

bias  is  increased  the  depletion  width  gets  larger,  which  in  turn 

increases  the  active  volume  of  the  device.  This  should  increase 

the  transient  current  due  to  the  depletion  component  by  a  factor 
1/3 

of  V  J ;  that  is , 

ip(t)  =  q  Acgwt  =  IQ  +  KV1/3  (162) 

where  I  is  the  zero  bias  depletion  component,  K  is  a  constant, 
and  V  is  the  magnitude  of  the  applied  voltage. 
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APPENDIX  III 


AFWL  600-KV  FLASH  X-RAY  FACILITY 

The  AFWL  flash  X-ray  system  consists  of  a  600-KV  Fexitron 
unit  manufactured  by  Field  Emission  Corporation.  The  character¬ 
istics  of  the  unit  are  summarized  in  table  XI.  A  double-walled 
rf  shield  room  houses  the  experimental  measuring  equipment.  By 
this  means  the  rf  noise  level  is  reduced  by  120  db.  The  X-ray 
beam  is  fired  through  a  special  window  in  the  shielded  room. 

This  window  is  shown  in  figure  37  •  The  beam  area  and  dose  rates 
are  shown  in  figure  38.  The  room  layout  is  shown  in  figure  5# 
Section  I.  The  energy  spectrum  of  the  system  is  shown  in  fig¬ 
ure  39  for  400-KV  X  rays.  A  graph  showing  the  conversion  of 
gamma  flux  equivalent  to  1  roentgen-hour  as  a  function  of  gamma 
energy  is  shown  in  figure  40. 

Table  XI 

Operating  Characteristics  of  the  600-KV 
Fexitron  Flash  X-ray  System 


Output  voltage  (peak) 

600  KV 

Output  current  at  peak 
output  voltage 

2000  amp 

Pulse  width  (at  half 
maximum  amplitude) 

0.2  nsec 

Pulse  rise  time 

50  nanosec 

Dose  rate  maximum 

«  2  x  10?  r/sec 
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Scale:  0.5  cm  =  1  inch 


m 


too  a\ 

H  tO^t 

•  I  • 

o  o  o 


I  I  I 


>  >  > 
XXX 


to 


Figure  38.  Collimated  Radiation  Beam  of  AFWL  600-KV  Flash  X-Ray  System 


Figure  59-  Energy  Spectrum  of  a 
400  kvp  Flash  X-Ray  Pulse 
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Gamma  Energy  (mev) 

Figure  40.  Gamma  Flux  Equivalent  to  1  Roentgen-Hour 
as  a  Function  of  Gamma  Energy 
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APPENDIX  IV 


THE  EFFECT  OF  TRANSIENT  RADIATION  ON  A 
GaAs  LIGHT-EMITTING  DIODE 

The  effect  of  transient  ionizing  radiation  on  a  General 
Electric  LED-10/  infrared,  light -emitting,  gallium  arsenide 
diode  was  experimentally  investigated.  The  infrared  light  from 
the  diode  was  monitored  using  a  newly  developed  EG&G  SD-100 
silicon  photodiode.  The  photodiode  has  a  spectral  response  from 
3500  to  11,300  angstroms  with  typical  rise  and  fall  times  of 
4  x  10~^  and  15  x  10~^  sec,  respectively.  The  typical  sensitiv¬ 
ity  of  the  photodiode  operating  at  9000  angstroms  was  0.25 
pia/iaw.  The  infrared  light -emitting  diode  had  a  spectrum  which 
peaked  at  approximately  8,970  angstroms  at  25°  C.  The  diode 
output  light  power  is  typically  0.3  mw  when  operated  at  a  con¬ 
tinuous  forward  current  of  100  ma. 

The  GE  diode  was  positioned  in  a  1-inch  diameter  collimated 
X-ray  beam.  It  was  connected  in  a  circuit  as  shown  in  figure  4l. 
A  51"°^™  resistor  and  a  maximum  coaxial  cable  length  of  twelve 
inches  were  used  to  minimize  the  viewing  circuit  time  constant. 

The  photodiode  was  aligned  with  the  output  beam  of  the 
light -emitting  diode  at  a  distance  of  approximately  0.375  inch. 
The  alignment  was  achieved  by  biasing  the  light -emitting  diode 
with  100  ma  DC.  Then  the  photodiode  was  moved  until  a  maximum 
voltage  was  250  microvolts. 

An  attempt  was  made  to  shield  the  photodiode  from  all  direct 
and  reflected  X  rays .  To  do  this,  a  hole ,  large  enough  to  accept 
a  TO- 18  size  transistor  can,  was  drilled  most  of  the  way  through 
a  2 -inch  thick  lead  brick.  The  remaining  0.1  inch  was  drilled 
to  the  same  size  as  the  lens  in  the  photodiode  can.  The  photo¬ 
diode  was  then  placed  in  the  hole.  The  brick  was  then  placed 
behind  another  2-inch  thick  lead  block.  The  end  result  was 
that  the  photodiode  was  shielded  from  direct  X  rays  by  5  inches 
of  lead,  and  from  reflected  X  rays  by  2  inches  of  lead.  See 
figure  42  for  the  physical  layout. 
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Coaxial  Cable 


Figure  4l.  Schematic  Diagram  of  Light  Emitting  Diode  Test  Circuit 


Figure  42.  Shielding  Layout  for  Light  Emitting  Diode  Tests 
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A  Polaroid  camera  attached  to  a  Tektronix  Type  551  oscillo¬ 
scope  was  used  to  monitor  the  transient  current  in  both  diodes. 
In  order  to  test  the  shielding,  a  test  was  performed  with  the 
GaAs  diode  removed.  The  resultant  picture  showed  a  low  noise 
level  and  no  apparent  radiation  response.  The  light-emitting 
diode  was  then  inserted  in  the  X-ray  beam  and  aligned. 

Tests  were  run  to  determine  the  effect  of  varying  bias 
voltage  on  the  response  to  transient  radiation.  As  the  forward 
bias  was  decreased,  the  transient  current  in  the  light-emitting 
diode  increased  as  shown  in  figure  43.  This  increased  transient 
current  with  decreasing  forward  bias  also  occurs  in  silicon 
diodes.  Therefore,  the  GaAs  diode  seems  to  respond  electrically 
to  ionizing  radiation  in  a  manner  similar  to  that  of  other 
diodes . 

The  increasing  transient  current  should  indicate  a  change 
in  light  emitted  by  the  diode;  however,  a  corresponding  increase 
was  not  detected  in  the  transient  current  of  the  photodiode. 

The  light-emitting  diode  was  then  moved  slightly  out  of  align¬ 
ment  with  the  X-ray  beam.  The  transient  response  of  the  photo¬ 
diode  remained  the  same  as  it  was  during  the  previous  tests. 

The  GaAs  diode  was  replaced  by  a  transistor  in  a  TO-18  can  and 
the  results  were  still  the  same  as  before.  Only  when  the 
transistor  was  removed  completely  was  the  transient  current  in 
the  photodiode  eliminated.  These  tests  indicated  that  the 
transient  current  in  the  photodiode  was  due  entirely  to  reflec¬ 
tions  from  the  packaging  material  that  surrounds  the  light- 
emitting  diode.  Tests  were  made  to  determine  the  effect  of 
varying  radiation  intensity  on  the  transient  current  in  the 
GaAs  diode.  Figure  44  shows  the  results  which  were  as  expected. 
The  peak  transient  current  increased  with  an  increase  of  radia¬ 
tion  intensity. 


All  Coaxial  Cabla 
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As  a  result  of  the  above  investigation  the  following 
conclusions  were  reached: 

(1)  The  GE  diode  is  relatively  insensitive  to  dose  rates 
of  less  than  10^  r/sec. 

(2)  A  more  sensitive  photodiode  or  higher  dose  rate  is 
needed  to  observe  any  possible  change  in  light  power  output. 

(3)  Operating  the  GE  diode  at  higher  forward  current  levels 
will  make  the  diode  more  radiation  hardened  to  ionizing  radiation. 

These  results  compared  favorably  to  those  obtained  by  other  in¬ 
vestigators  (references  18  and  19)* 
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